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aligned  with  the  total  velocity  vector,  the  z axis  in  the 
aircraft  plane  of  synmetry  pointing  down.-  and  the  y aids 
orthogonal  to  the  x and  z axes  and  positive  out  the  right 
wing 

Zi  l^(  )/^i]  /■ 

Controlled  element  transfer  function 

Pilot  describing  function  as  particularized  by  subscript 
Angle  of  attack 
Sideslip  angle 
Fli^t  path  angle 

Surface  deflection  particularized  by  subscript 
Transfer  function  denominator 
Damping  ratio  particularized  by  subscript 
ConventiOTial  perturbed  attitude  angles 

Conventional  Euler  angles  between  gravity- oriented  inertial 
axis  and  aircraft  body  axis 

Angle  between  thrust  axis  and  edrcraft  centerline  x axis 

Real  part  of  Laplace  transform  variable 

Nat’ural  frequency  of  the  dencmiiAtor  or  nimerator  root 
particularized  by  the  siibscript 
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SECTIGir  1 


UnEODUCTIQN 


Modern  and  advanced  military  aircraft  often  fly  at  high,  near-stall 
angles  of  attack  in  tactical  operations  and  "normal"  maneuvering.  In  fact, 
superior  flight  characteristics  in  this  marginjJ.  regime  are  key  effective- 
ness factors  in  the  world  of  fighter  and  fighter-bomber  cai5)etitions  and 
contests.  They  are  also  central  issues  in  survival,  since  slight  miscues 
at  these  boundary  conditions  can  lead  to  out- of- control  departures  and  spins. 
High  performance  aircraft  are  often  very  difficult  to  recover  from  such 
situations,  so  the  loss  of  control  can  be  very  dangerous. 

Unfortunately,  understanding  of  high-angle-of-attack  flight  is  not 
cotipatible  with  its  level  of  operationsd  importance.  A variety  of  uncon- 
trolled hi^-angle-of-attack  dynamic  phenomena  have  been  observed  and  cata- 
logued; some  have  been  analyzed  for  underlying  physical  causes.  However, 
there  is  little  present  appreciation  and  quantification  of  the  impact  of 
root  deficiencies  in  high-angle-of-attack  aircraft  dynamics  on  pilot- 
controlled  flying  qualities  just  before,  dur..ng,  and  just  after  stall,' 

departure.  For  instance,  analyticsd  methods  for  efficient  treatment  of 
the  total  pilot /aircraft  syc'  em  have  not  been  applied  specificsilly  to  high- 
angle-of-attack  flying  qualities  prjblems.  Most  criteria  have  evolved  from 
open-loop  (airframe  sQ-one)  considerations  yet  departure  itself  represents  a 
dividing  line  in  the  nature  of  the  approaches  necessary  to  modeling  and 
ansJ.ysis  requirements.  Before  departure  the  control  situation  is  primarily 
regulatory,  closed  loopj  after  departure,  the  control  situation  through 
recovery  is  ideally  time-optimal  programmed  (i.e.,  largely  open-loop).  In 
view  of  the  newly  re-recognized  importance  of  high-angle-of-attack  flight, 
and  because  the  application  of  the  genereil  pilot/vehicle  systems  analysis 
methods  to  the  attendant  closed-loop  control  problems  is  very  promising, 
the  time  for  a breakthrough  in  understanding  was  considered  to  be  at  hand. 


The  specific  goals  of  this  program  were: 

• Enhanced  understanding  of  hi^  angle- of- attack  ppen- 
and  closed-loop  dynamic  phenomena 

• Identification  of  key  effective  controlled  element 
flying  qiiality  parameters 

• Improved  appreciation  for  interactions  between  hig^i 
angle- of- attack  flying  qualities,  controlled  element 
dynamics,  and  pilot  technique 

• definition  of  multi- input,  multi-output  pilot  models 
; lippropriate  in  the  stall/departure/rec every  phases 

' of  flight 

• Develqpnent  and  validation  of  quantitative  design 
guides  and  criteria  relevant  to  stall/departure/ 
recovery 

• Formulation  of  an  automated  closed-loop  pilot/vehicle 
digital  design  analysis  methodology 

• Recfmmendations  for  further  work  and  approaches  to 
be  ;aken  by  the  AFFDL 

With  the  exception  of  the  automated  digital  closed-loop  pilot/vehicle 
design  analysis  method,  all  of  these  goals  have  been  attained. 

Attention  has  been  focused  on  dqjarture  (e.g.,  nose  slice,  wing  rock) 
and  recovery  as  opposed  to  spin  entry  and  recovery.  The  approach  has  been 
to  analytically  identify  and  validate  key  flying  quality  parameters,  to 
qualitatively  assess  and  verify  these  throu^  piloted  fixed-base  simulations, 
and  to  measure  and  quantify  effects  on  key  dynamic  parameters  of  pilot 
describing  functions. 


Flying  qualities  are  influenced  by  open-  and  closed-loop  psurameters. 

Most  past  studies  (e.g..  Refs.  1-4)  have  concentrated  on  the  open-loop 
lateral-directional  parameter  contributions  (e.g.,  Np,  Xp)  sdthough 

some  recognition  has  been  given  to  the  closed-loop  parameters  and 
(e.g..  Refs.  5 and  6).  These  efforts  have  also  focused  on  symmetric  (p  = 0) 
aircraft  trim  conditions. 
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[ The  key  parameters  investigated  in  this  program  are  identified  in 

* Table  1,  Since  the  pilot  is  generally  in  the  loop  when  departure  occurs, 

we  place  most  emphasis  on  key  closed-loop  parameters  and  then  relate  these 
to  open-loop  parameters  (which  in  this  case  are  aircraft  stability  derivatives). 
( Several  of  these  parameters  are  new  and  result  primarily  from  recognizing  that, 

j at  high  a,  it  is  conmon  for  sidelip  to  exist  either  intentionally  (e.g.,  rudder 

maneuvering)  or  unintentionally  (e.g.,  adverse  aileron  yaw,  mistrim,  etc.).  As 
a direct  result  the  new  open-  and  closed-loop  parameters  have  been  identified 
! which  evolve  from  longitudiml-lateral  coupling  associated  with  unsymmetric 

fli^t,  Uiese  new  parameters  relate  to  the  zeros  (i/Tq^,  ojg)  of  the  pitch 
attitude  transfer  function  reaching  unfavorable  locations  at  angles -of- at tack 

I below  that  for  stall  and  at  relatively  small  sideslip  angles.  Conventional 

feedback  of  pitch  attitude  (or  rate)  to  elevator  (either  automatic  or  manual) 
then  produces  a coupled  longitudinal- lateral  divergence  known  as  nose  slice. 

S! 

r 

I 'Hie  Section  II  analysis  is  focused  primarily  on  the  A-7  aircraft  which 

■ exhibits  a severe  nose  slice  departure  characterized  by  large  rapid  yaw 

followed  by  rapid  roll.  Limited  analysis  and  simulation  of  the  F-U  wing 
rock  is  also  reported.  Physical  insight  to  the  nose  slice  phenomena  is 
I achieved  by  tracing  through  the  equation  of  motion,  effective  stability 

’ derivative,  transfer  function,  and  root  aerodynamic  flew  relationships. 

► Eighth  order  transfer  functions  are  obtained  for  all  state  variables  for  a 

? sweep  of  a and  P in  the  vicinity  of  stsll.  It  is  found  that  the  pitch 

TAHLE  1 

KEY  PARAMETERS  INVESTIGATED 


DEPARTURE 

CHARACTERISTIC 

PARAMETERS  INVESTIGATED 

CLOSED-LOOP 

OPEN-LOOP 

Nose  Slice 

1/^05  , 0)0 

o)SR  ^ 

Np 

*^dyn 

Wing  Rock 

?d“d  ^ ‘^p/^d 

nonlinear  Np 
Np  , Xp 
Mp 

attitude  numerator  undergoes  a marked  change  vioh  non- zero  p In  vhlch  one 
root  becomes  non-minimum  phase  (moves  into  the  right  half  plane  of  a con- 
ventional root  locus  plot).  IMs  portends  a closed  loop  instability.  The 
cause  of  this  right  half  plane  zero  is  traced  through  the  9*9  equation 
of  motion  matrix  to  the  dominant  coupling  terms.  Two  of  these,  (jCi  ahd  N^, 
are  aerodynamic  idiile  the  other  two,  Zp  and  Zy,  are  kinematic  coupling. 

The  potential  closed-loop  unstable  mode  is  identified  as  a first  order 
lateral-directional  divergence  despite  its  appearance  in  a longitudinal 
numerator.  The  phenomenon  is  verified  by  open-  and  closed- loop  time 
histories  obtained  from  nonlinear  six  degree-of-freedom  digital  and  five 
degree-of- freedom  analog  simulations  using  nonlinear  aerodynamic  coefficents. 
The  re stilts  are  further  validated  in  a fixed  base  piloted  simulation  using 
the  five  degree-of-freedom  nonlinear  airframe  model.  Section  III  contains 
a description  of  the  simulation  including  head-up  display,  tracking  tasks 
and  pilot  parameter  measurements. 

Nine  aircraft  dynamic  configurations  are  enployed  in  the  departure/ 
recovery  flying  quality  assessment  simulation.  These  and  the  simulation 
results  are  analyzed  in  Section  IV,  An  in-trail  target- tracking  task  was 
flown  by  two  USAF  test  pilots  experienced  in  the  A-7,  This  task  led  the 
subjects  into  a level,  1 g,  stall  departure.  Initial  values  of  each  of 
the  parameters  of  Table  1 were  varied  to  provide  the  nine  departure/recovery 
cheuracteristics  ranging  from  mild  to  severe  nose  slice  and  with  one  provid- 
ing a tendency  to  wing  rock. 

Analytic  predictions  are  validated;  however,  resolution  of  pilot  commentary 
required  additional  analysis  which  led  to  further  identification  and  separation 
of  key  parums+er  effects  and  possible  analytic  predictive  techniques.  Departure 
characteristics  are  found  to  be  determined  by  the  closed-loop  parameters  of 
Table  1.  Recovery  techniques  are  found  to  be  in  agreement  with  A-7  flight 
results  and  the  preferred  method  is  release  of  all  controls;  however,  in  some 
cases  application  of  forward  stick  is  acceptable.  Therefore  recovery  charac- 
teristics are  strongly  influenced  by  the  open- loop  parameters  of  Table  1 . 
Departure  characteristics  and  piloting  technique  also  highly  influence  recovery. 
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Pilot/vehicle  describing  function  measures  and  pilot  dynamic  models 
obtained  for  three  incipient  departure  conditions  are  summarized  in  Section  V. 
The  task  is  again  representative  of  a situation  in  which  aircraft  stall  is 
inadvertently  approached  and  the  pilot  is  attempting  to  maintain  track  without 
inducing  aircraft  depeirture.  Describing  function  measures  were  obtained  for 
raultiloop,  pitch,  roll,  and  heading  tracking  tasks.  The  crossover  model  was 
found  to  be  valid  for  each  loop  closure.  The  precision  pilot  model  provided 
an  excellent  fit  to  the  data  except  that  in  pitch  the  crossover  was  close  to 
the  longitudinal  short-period  mode  and  required  the  pilot  to  adopt  a second- 
order  lead  equalization  approximately  the  inverse  of  the  longitudinal  short 
period.  The  gain  and  low-frequency  lag- lead  equalization  employed  in  the 
roll  loop  closure  was  found  to  differ  depending  upon  vrtiether  the  roll  loop 
was  an  outer  (roll  tracking)  or  inner  (heading  tacking)  loop.  Coupling 
between  the  airframe  lateral  and  longitudinal  dynamics  was  found  to  increase 
the  pilot  workload  but  did  not  appreciably  affect  the  equalization  enployed. 

The  one  goal  not  achieved  — the  autoinated  digital  closed-loop  pilot/ 
vehicle  design  analysis  method  — was  set  forth  as  a continuation  and  multi- 
input, multi- output  extension  of  the  vsu-ious  "paper  pilot"  models  (Refs.  7, 

8,  and  9)*  As  such  it  was  anticipated  to  draw  on  the  existing  data  base 
relating  pilot  ratings  with  first-order  lead  generation  requirements  and 
task  performance  measures  in  establishing  a rating  functional.  However,  this 
preconceived  concept  did  not  fully  recognize  the  grcssness  of  the  maneuvers, 
pilot  control  activity,  and  performance  measures  associated  with  stall/ 
departure/recovery  phases  of  flight.  DQ>arture  is,  by  definition,  the  attain- 
ment of  aircraft  motion  uncomraanded  by  the  pilot,  i.e.,  loss  of  control.  The 
motion  is  generally  large  and  rapid.  It  is  found  in  this  study  that  departure 
(nose  slice  or  wing  rock)  is  severely  aggravated,  if  not  caused,  by  closed- 
loop  pilot- vehicle  interaction.  With  the  possible  exception  of  holding 
forward  stick,  recovery  is  best  achieved  by  the  pilot  eleasing  all  controls. 
Any  other  control  activity  during  recovery  greatly  aggravated  the  departure/ 
recovery  maneuver.  The  pilot  ratings  received  for  closed- loop  recovery  were 
therefore  much  higher  (worse)  than  for  open-loop  airframe-alone  recovery. 

Thus  a pilot -centered  function  is  not  applicable  to  recovery. 
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In  the  pre-departure  phase  \rtiere  closed-loop  considerations  aure  valid, 
it  was  found  the  the  pilot  adopts  a seccmd-order  lead  in  pitch  attitude 
control  — apparently  in  an  atten5>t  to  cancel  short  period  pitch  attitude 
excursions  which  might  trigger  departure.  Tliis  second-order  ccn^jlex  lead 
generation  has  been  observed  before  (Refs,  10  and  11);  however,  there  are 
insufficient  pilot  rating  data  to  establish  a functior.fJ..  The  only  other 
known  rating  data  is  in  Ref,  10  and  it,  together  with  the  data  obtained 
herein,  provides  essentially  no  rating  spread.  The  vast  amount  of  data 
obtained  for  first-order  lead  generation  is  net  applicable.  Thus  at  this 
point  it  appears  a paper  pilot  design  function  is  somewhat  premature. 

Nevertheless,  a pilot  model  of  a different  sort  has  been  obtained  for 
the  high-*angle-of-attack,  pre-departure  flight  phase.  However,  much  more 
data  is  required  before  trends  can  be  'stablished  between  pilot  parameters, 
pilot  ratings,  and  performance  measiires. 

A key  problem  with  ratings  exists  because  loss  of  control  is  involved. 
That  is,  accomplishment  of  the  primary  task,  whatever  it  may  be,  has  resulted 
in  the  aircraft  approaching  stall/departure.  If  the  task  is  pursued,  it 
presumably  will  result  in  departure  — loss  of  control.  Departure  prevention 
then  requires  abandonment  of  the  primary  task  and  shifting  full  priority  to 
safety  of  flight.  Thus  Level  5 flying  qualities*  are  involved  and  the 
problem  arises  of  establishing  the  Level  ^ boundary,  i.e.,  loss  of  con- 
trol. On  the  Cooper-Harper  ^'Ref.  15)  scale  this  is  the  boundary  between 
9 and  10  (there  is  no  9.5).  If  an  airplane  will  depart,  then  there  is  no 
question  regarding  the  existence  of  the  Cooper-Harper  boundary,  but  merely 
a question  of  hew  much  of  the  potential  fli^t  envelope  must  be  traded  for 
flight  safety.  Thus  a mission  performance,  usable  load  factor,  flight  safety, 
or  seme  related  index  becomes  paramount. 


•From  Ref.  12,  Level  5: 

"Flying  qualities  such  that  the  airplane  can  be  controlled 
safely,  but  pilot  workload  is  excessive  or  mission  effective- 
ness is  inadequate,  or  both.  Category  A Flight  Phases  can  be 
terminated  safely,  and  Category  B and  C Flight  Phases  can  be 
completed. 


•V 


In  summELiy,  the  study  reported  herein  was  intended  to  unify  into  a 
systematic  framework  a wide  variety  of  empirical,  analytical,  and  theo- 
retical data  into  analytical  models  which  describe  high-angle-of-attack 
aircraft  dynamics  and  flying  qualities.  Besides  providing  the  basis  for 
enhanced  'mowledge,  this  work  also  was  intended  to  develop  and  validate 
quantitative  design  guides  and  associated  flying  qualities  and  fli/^t  test 
criteria  relevant  to  the  stall/departure/recovery  regime.  Spin  aspects  are 
specifically  excluded.  However,  because  of  the  ccnplexity  involved  in  non- 
linear aerodynamics  and  multi-input,  multi-output  pilot  characteristics, 
the  effort  has  been  a con^jrehensive  and  systematic  first  cut.  It  is  antici- 
I>ated  that  the  insist  and  initial  flying  qualities  results  provided  will 
serve  in  planning  a longer-term  program  which  attacks  the  overall  problem 
of  sta?J.s  and  spins. 
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Siamia  AIULYSIS  HJEHTIFZCAIIDN  of  KET  SARANEIBtS 
A.  slice 

This  section  presents  the  analytical  development  of  and  possible 
physical  explanation  for  new  open-  and  closed-loop  parameters  which  appear 
to  underlie  the  nose  slice  departure  of  the  A-7  aircraft.  The  initial 
analysis  is  based  on  a nonlinear  six-degree-of-freedom  digital  model  which 
incorporates  lookup  table  aerodynamic  coefficients  as  functions  of  a,  p, 
and  6e«  A succession  of  fixed  operating  point,  partial  derivative 
expressions  and  linearized  transfer  function  evaluations  over  a range  of 
angles  of  attack  and  sideslip  provided  the  initial  clues  to  potential  key 
parameters.  Open-  and  closed- loop  time  responses  were  then  obtained  for  a 
nonlinear  six-degree-of-freedom  digital  and  a simplified  five-degree-of- 
freedom  analog  simulation.  These  models  are  validated  via  comparisons  with 
actxial  fli^t  traces  of  an  A-7  departure.  The  closed-loop  nose  slice  diver- 
gence phenomena  is  traced  throu^  effective  stability  derivatives  and  related 
to  previous  aerodynamic  flow  studies. 

1 . Trim  Conditions 

Low  Mach  aero  data  developed  by  the  manufacturer  were  supplied  by  the 
USAF  Flight  Dynamics  Laboratory.  These  data  cover  the  range  0 < a < 90° 

-20°  < P < +20°.  The  data  supplied  by  the  '^F  was  in  stability  axes,  suid 
was  trsuisforraed  into  body  centerline  axes  fo’-  this  analysis  emd  simulation. 
Volume  II  contains  the  body  axis  data. 

Trim  conditions  for  this  auialysis  and  simuiation  arj  as  follows 

W = 22,299  lb  = 0 

h = 15,200  ft  Landing  ge&r  — Up 

*cg  = 29.6^  MAC  Haps  — Up 
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Tlie  velocity  was  varied  to  obtain  an  angle-of-attac^  range  10  < a < 2h  deg. 
Since  the  clean  aircraft  stalls  at  a = 21  deg,  this  range  provided  an  an5>le 
view  of  the  vehicle  dynamic  characteristics  before,  through,  and  following 
stall.  It  was  subsequently  determined  that  an  a range  of  from  17  to  24  deg 
is  adequate  to  demonstrate  departure /recovery  characteristics  for  this  aircraft 

The  linearized,  fixed  operating  point  derivatives  and  transfer  function 
factors  wre  obtained  throu^  use  of  standard  digital  programs.  These  pro- 
grams utilize  the  equations  of  motion  summarized  in  Fig.  1 and  the  nonlinear 
lookup  table  aero  data  of  Appendix  II.  The  programs  trim  the  aircraft  to  the 
specified  flight  condition,  perturb  each  of  the  state  variables  to  obtain  the 
partial  (stability)  derivatives,  and  compute  the  specified  transfer  function 
factors . 

The  initial  analysis  at  p = 0 showed  rather  conventional  uncoupled  dynamic 

characteristics  and  provided  no  clues  to  a nose  slice  type  departure.  The 

tabulated  data  in  Appendix  III  show  that  dutch  roll  frequency  decreases  and 
damping  increases  as  angle  of  attack  is  increased.  The  conventional  roll 
subsidence  and  spiral  modes  couple  into  a lateral-oscillation  (hereafter 
termed  "lateral  phugoid")  that  eventually  becomes  diveigent  between  20  and 
21  deg  angle  of  attack.  Longitudinal  dynamics  are  well  behaved  with  no  indi- 
cation of  pitch-up  or  other  longitudinal  departure  characteristics. 

Because  of  these  "negative"  results  for  p = 0 and  because  previous  exposure 

to  the  A-7  (Ref.  14)  iniicutei’  that  some  moiels  ;ia1  a pronounced  tendency  to 

change  in  directional  trim  during  accelerating  or  decelerating  flight,  the 
effect  of  non- symmetrical  flight  was  investigated.  Six-degree-of-freedom 
transfer  function  factors  for  "hree  sideslip  values  at  a = 19  deg  are  shown 
in  Table  2.  It  may  be  noted  from  the  characteristic  equation  (A)  factors  that 
the  open-10'7  longitudinal  modes  are  little  affected  by  sideslip.  The  major 
influence  is  on  the  open- loop  lateral  directional  modes.  Increasing  sideslip 
causes  the  roll  subsidence  and  spiral,  to  couple  into  a lateral  oscillation 
much  as  does  increasing  angle  of  attack.  The  damping  of  the  lateral  phugoid 
decreases  with  increasing  P and  goes  unstable  between  6 a.id  15  deg  of  side- 
slip. Comparing  the  ducch  roll  and  lateral  phugoid,  it  appears  that  there  is 
an  interchange  in  damping  between  these  modes.  That  is,  as  sideslip  increases, 
damping  of  the  dutch  r^j.1  increases  while  that  of  the  lateral-phugoid  decreases 
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Figure  i . Equation  Summary 


TABLE  2.  SURVEY  OF  TRANSFER  FUNCTIONS  WITH  CHANGE  IN  p 


On  the  other  hand.  It  may  be  noted  from  the  lateral-directional  numerator 
factors  that  sideslip  has  little  effect  on  the  zeros  of  the  lateral  or  longi- 
tudinal degrees  of  freedom.  However,  there  is  a significant  difference  for 
the  pitch  rate  and  pitch  attitude  to  elevator  numerators.  These  show  l/Tei 
to  be  small  but  negative  (non-minimum  i^se)  at  zero  sideslip  and  to  rapidly 
move  to  much  larger  negative  values  as  sideslip  increases.  At  the  same  time 
1 /T02,  which  is  initially  positive,  moves  to  much  larger  positive  values 
with  increasing  sideslip.  In  fact,  the  movement  of  l/Te^  and  I/T02  are  pre- 
cisely the  opposite  to  that  of  the  denominator  roots  l/Ts  and  1/Tp.  The  pitch 
numerator  roots  move  away  from  each  other  as  sideslip  is  Increased  whereas 
the  denominator  roots  move  towards  each  other  and  coalesce  into  the  lateral- 
phugoid  mode.  The  large  negative  value  of  l/To^  indicates  a potential  insta- 
bility vqpon  closure  of  either  the  pitch-rate-  or  pitch-attl tude-to-elevator 
loops.  This  is  of  considerable  concern  since  we  would  assume  the  pitch  attitude 
loop  would  be  tightly  closed  by  the  pilot  at  hi^  angle  of  attack.  Further- 
more, this  would  seem  to  indicate  a longitudinal  departure  mode  Instead  of  a 
lateral  departure  mode.  Attention  was  therefore  turned  to  the  matrix  formu- 
lation to  identify  key  off-diagonal,  or  coupling,  terms  which  might  be  influenc- 
ing these  shifts  in  numerator  and  denominator  roots. 

2.  Identification  of  Key  BarasKters 

Figure  2 presents  the  nine  by  nine  matrix  (three  body  axis  moments,  three 
flight  i>ath  forces,  and  three  Suler  angle  transformation  equations)  for 
coupled,  non- symmetric  flight  obtained  from  the  partial  derivative  expansion 
of  the  Fig.  1 equations.  Only  off-die^sonal  terms  which  analysis  determined 
to  have  a significant  effect  on  the  transfer  function  are  identified.  For 
the  A-7,  the  major  coupling  is  provided  by  the  terns  within  the  heavy  border- 
lines. Two  of  these,  X<x  and  N^,  are  aerodynamic  and  two,  Pq 
Po  sin  oo,  are  nonlinear  kinematic  terms.*  [ihe  coupling  derivative  Mp  is 

« 

Tile  earliest  known  discussion  of  this  coupling  phenomena  was  documented 
in  1916  (Ref.  15)  for  a six-degree-of- freedom  analysis  of  the  B.E.-2  biplane 
in  which  the  octals  were  factored  by  hand:  This  aircraft  exhibited =Xp; 

and  ^^  = Hie  conclusion; 

"It  would  be  difficult  to  build  a machine  in  other  respects  satis- 
factory which  would  remain  con^letely  stable  under  very  rapid 
side-slipping.  A small  change  either  in  Increase  or  decrease  of 
rudder  area  would  be  liable  to  cause  either  a spiral  dive  or  an 
(continued  on  following  pe^e) 
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essentially  zero  for  the  A-7  as  computed  from  the  available  data.]  Also 
shown  in  each  matrix  element  is  the  magnitude  of  the  term  evaluated  at 
do  = l8.8  deg  'and  Po  = 6 deg. 

Note  in  the  Z equation  that  the  effective  derivative  Zp  = cos  is 
roughly  one-third  the  value  of  while  Zj.  = Pq  sin  oq  is  about  one-tenth 
of  Z^.  Thus,  one  might  expect  these  coupling  terms  to  alter  the  pitch 
numerator  zero  1 /T02  since  for  the  uncoupled  case  1 /T0g  = -Zw.  In  the  roll 
equation, Xa  is  about  three-quarters  the  magnitude  of  Xp  and  of  opposite 
sign;  thus,  for  this  trim  condition  (positive  p),  an  increase  in  a results 
in  a moment  which  tends  to  cancel  the  roll  static  stability.  In  the  yaw 
equation  is  about  eig^t  times  the  magnitude  of  N^  and  of  the  same  sign. 

Both  ni,  and  Np  are  negative  and  would  indicate  a strong  directional  diver- 
gence characteristic  with  increasing  a or  p. 

The  effect  of  the  off-diagonal  terms  on  pitch  attitude  transfer  function 
pole-zero  locations  is  demonstrated  in  Fig.  5-  Figure  3a  shows  a completely 
uncoupled  six-degree-of- freedom  case,  for  reference,  in  which  the  two  lateral- 
directional  moles  (cD^  and  ojgg)  have  cancelling  pole-zero  dipoles  as  would  be 
expected.  The  longitudinal  poles  and  zeros  reflect  the  near  steLLl  condition 
in  which  1 /Teg  becomes  small  and  negative  1/T0i  might  be  expected.  Figure  3d 
presents  the  pole-zero  locations  for  the  completely  coupled  six-degree-of- 
freedom  case  which  shows  the  poles  to  be  little  affected  by  coupling,  whereas 
a major  shift  occurs  in  zero  locations.  The  most  significant  movement  is  in 
the  two  real  zeros  because,  as  will  be  shown  later,  closure  of  the  0 -»>  6e 
loop  will  drive  one  closed- loop  root  toward  the  right  half  plane  (RHP)  zero 
and  system  instability.  Of  secondary  concern  is  the  separation  of  the  dutch 

(continued  from  preceding  page) 

instability  of  spin,  necessarily  followed  by  a nose  dive  either 
with  or  without  previous  stalling,  in  order  to  obtain  the  neces- 
sary speed  to  overcome  the  instability.  The  large  decrease  in 
the  damping  factor  of  the  phugoid  oscillation  would  tend  to 
accelerate  this  performance." 

A similar  development  to  the  above  but  eti5)loying  a five-degree-of-freedom 
model  is  presented  in  Ref.  16  for  a swept-wing  transport  aircraft.  The  same 
coupling  terms  were  found  to  dominate  at  much  lower  angle  of  attack  but  com- 
parable sideslip  range  (e.g.,  6 to  i;:;  deg).  In  the  Ref.  l6  study,  the  coupling 
primarily  produced  an  oscillatory  aerodynamics /kinematic  interchange  in  short- 
period  energy  between  p and  a via  La- 
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roll  dipole  pair  which  indicates  considerable  modal  response  excitation  via 
elevator  (i.e.,  longitudinal  to  lateral  coupling). 

The  influence  of  the  individual  pair  of  coupling  terms  is  identified  in 
Figs.  3b  and  3c.  In  Fig.  3b  the  Z equation  off-diagonal  terms  are  set  to 
zero^  in  Fig.  3c  the  L and  N equation  off-diagonal  terms  are  removed.  Both 
result  in  similar  influences  on  the  various  zeros  and  indicate  these  "effective” 
derivatives  must  occur  in  combined  or  multiplicative  form  in  the  transfer 
function  numerator. 

As  an  aid  in  identification  of  the  modes  reflected  by  the  poles  emd  zeros 
of  Fig.  3,  the  longitudinal  dynamics  were  reduced  to  a short-period  approxi- 
mation by  deletion  of  the  X equation,  which  usually  eliminates  the  I/Tq^  zero 
and  replaces  the  conqplex  jtogoid  pole  with  a first-order  pole  at  the  origin. 
However,  Fig.  3e  shows  the  zeros  previously  identified  as  l/Tei  and  1/T02  to 
remain  unchanged  from  the  complete  six-degree -of- freedom  case.  The  complex 
zero  previously  identified  as  a>sg  has  become  a first-order  zero  near  the 
origin.  The  phugoid  mode  is  transformed  into  a first-order  pole  at  the  origin 
as  expected.  Because  the  pole-zero  configuration  of  Fig.  3e  reflects  coupled 
lateral-longitudinal  modes  (compare  with  Fig.  3d),  the  real  zero  in  the  right 
half  plane  previously  labeled  I/T91  will  be  identified  as  l/Te^  in  subsequent 
analysis  since  this  is  a new,  coupled  lateral-longitudinal  mode. 

In  a further  attempt  to  identify  the  pitch  attitude  zeros,  the  matrix  for 
the  five-degree-of-freedom  pitch  attitude  numerator  was  expanded  in  litereJ. 
terms  and  the  polynomial  coefficients  obtained.  Approximations  containing 
only  the  most  significant  terms  are  presented  in  Eq.  1 . 

N6g  = MBg[As5  + Bs*^  + Cs^  + Ds^  + Es  + F]  (l ) 

where  (coupling  terms  are  underlined): 

A = 1 

B = z;  + Xp  + (n;  + Yv) 

C = X^(Zw  + N^)  sin  a 

D = “Xp[(g/tJo)  cos  9 — (Zy  + Np)  sin  a]  ~Np(Zy +<p)  cos  a 
E = Xp(g/Uo)  cos  9[(Zw4N;)-(Z6e/tfee)Mw]+2p(^(^g-N^) 

F = ^g(g/Uo)  cos  9N;Zw 


TR- 1033-1 
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It  may  be  observed  that  the  off-diagonal  coupling  terms  are  multiplicative 
and  primarily  influence  the  C and  E coefficients.  An  attempt  was  made  to 
factor  this  equation  into  literal  expressions  for  the  iraots.  However,  this 
task  proved  too  time  consuming  and  was  abandoned. 

a.  Single-Loop  Control 

A single-loop  system-survey  for  elevator  control  of  pitch  attitude  with 
the  six-degree-of-freedom  coupled  airframe  in  non -symmetrical  flight  is  shown 
in  Fig.  4.  The  transfer  function  is  shown  in  the  upper  left.*  The  root  locus 
in  the  top  right  of  the  figure  reflects  root  migrations  for  a pure  gain  closure. 
Note  that  the  roots  starting  at  u)gj^  rapidly  move  to  the  real  axis  and  then 
split  into  two  real  rootsj  one  of  which  moves  towards  1 /T02,  the  other  moves 
towards  l/Te^.  The  rapidity  of  the  movement  of  these  closed- loop  poles  towards 
the  zeros  is  demonstrated  by  a 3ode-siggy  plot  in  the  bottom  half  of  Fig.  U. 

The  heavy  solid  and  dashed  lines  of  the  Bode  correspond  to  the  path  of  the 
closed-loop  roots  along  the  real  (o)  axis  in  the  root  locus  above.  As  the 
loop  gain  is  increased,  the  cortplex  poles  emanating  from  cugj  meet  the  real 
axis  at  the  apex  of  the  solid  curve  in  the  Bode-siggy  plot.  Further  increase 
in  gjain  moves  one  closed- loop  root  to  a lower  frequency  or  towards  the  origin 
while  the  other  root  moves  to  higher  frequency  and,  at  very  high  gain,  asymp- 
totically approaches  the  1 /T02  zero  at  0.866  rad/sec.  The  root  that  goes  toward 
the  origin  passes  into  the  right  half  pleuie  as  shown  in  the  root  locus.  This 
is  represented  in  the  Bode-siggy  bj  the  dashed  line  which  reflects  the  mirror 
image  of  the  closed-loop  pole  asyraptoticedly  approaching  the  l/T0j  zero  at 
-0.5  rad.  If  the  pilot  is  to  achieve  effective  control  of  pitch  attitude,  he 
must  close  the  loop  so  the  gain  line  lies  below  the  low-frequency  asymptote 
of  the  Bode  plot.  It  is  obvious  that  this  then  results  in  a closed-loop  pole 
in  the  right  half  plane.  If  the  pilot  closes  the  loop  so  that  unity  gain 
"crossover"  is  achieved  in  the  region  of  1-5  rad/sec,  which  covers  the  range 
of  usual  "loose"  to  "tight"  piloted  pitch  attitude  control,  it  may  be  seen 
that  the  closed-loop  poles  will  lie  very  close  to  the  open- loop  zeros.  For 
example,  a unity  dc  gain  provides  a crossover  between  1.5  and  2.5  rad/sec  and 


“Throughout  tuis  report  a shorthand  notation  is  used  in  presenting  transfer 
functions:  (s  + a)  = (a)j  [s^  + 2^cds  + u^]  = [^,  cu]. 
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closed-loop  roots  at  -0.58  and  -K).66  rad/sec.  Ihe  resulting  first-order 
divergence  has  a time  constant  of  about  3*6  sec. 

The  difference  in  open-  and  closed-loop  responses  may  be  observed  from 
time  traces  of  the  six-degree-of-freedom  digital  simulation  using  nonlinear 
aerodynamic  lookup  tables.  Figure  3 shows  the  open- loop  lateral  and  longi- 
tudinal time  traces  for  the  A-7  trimmed  at  ein  initial  17.3  deg  a and  0 deg  3 
and  a ramp  6e  input  of  1.2  deg/sec.  An  aileron  doublet  is  sdso  introduced 
at  4 sec  to  provide  some  excitation  of  the  latered  modes.  The  traces  indicate 
that  a slow  directional  divergence  starts  at  approximately  8 sec  and  2k  deg  a. 
It  then  couples  with  p (sind  a)  to  form  a divergent  oscillation  with  a 
period  of  approximately  11  sec.  This  will  be  shown  in  subsequent  sections  to 
I be  an  unstable  lateral-phugoid  mode.  If  the  aircraft  were  trimmed  at  a small 

initial  p the  open-loop  response  would  remain  nearly  the  same  but  the  diver- 
gence would  develop  at  a sli^tly  lower  a because  of  the  combined  effect  of 
P and  a on  the  lateral  phugoid  as  noted  previously. 

Figure  6 shows  the  system  response  when  the  pitch  attitude  loop  is  closed 
I with  unity  gain  (as  in  Fig.  U)  and  a step  0c  of  0.01  rad  is  introduced.  Ihe 

i aircraft  is  initiadly  tri-med  for  steady  fli^t  at  a = 18.8  deg,  p = 6 deg,  and 

f 

■ cp  = 5 deg.  The  predicted  first-order  divergence  is  seen  to  start  immediately 

I and  to  dominate  the  r,  \lr,  'p,  9,  and  a traces.  The  divergence  time  constant  is 

' approximately  3*6  sec,  as  measured  in  the  r,  7,  0,  and  a traces,  consistent 

t with  the  closed- loop  right  half  pleuie  root  location  predicted  in  Fig.  k.  By 

the  time  2k  deg  a is  reached  the  aircraft  is  well  departed.  These  responses 
verify  the  closed-loop  mode  to  be  coupled  longitudinal-latereLL  motion.  They 
also  appear  to  match  the  description  of  the  classic  A-7  nose  slice  departure 
described  in  Ref.  17  as  an  initial  hi^  yaw  rate  followed  by  a hi^  roll  rate. 
Ihe  traces  of  Fig.  6 show  that  at  9 sec,  which  is  just  before  the  hij^  roll 
rate  (departure)  develops; 

r = — 15  deg/sec  (extrapolated) 
p = —*+.9  deg/sec 
P = +13.3  deg 
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■Riey  also  show  that  the  initial  longitudinal  comnand  excites  considerable 
roll  rate.  It  will  be  shown  later  that  both  the  longitudinal  short  period 
and  the  lateral  dutch  roll  contain  considerable  roll  rate  response.  These 
two  modes  then  result  in  the  ’’beat  frequency"  shape  of  the  roll  rate  time 
trace. 

b.  Multiloop  Control 

Lateral  control  is  used  very  sparingly  at  hi^  angles  of  attack  because 
it  tends  to  aggravate  directional  departvure.  Nevertheless,  it  is  pertinent 
to  assume  that  a pilot  would  tend  to  close  a simple  (pure  gain)  bcuik-angle- 
to-aileron  loop,  shown  in  Fig.  7,  under  the  lateral-directional  departure 


circumstances  indicated  in  Fig.  6.  The  use  of  a second  control  stirface  li^iut 
(e.g.,  roll  control  via  aileron)  with  the  coupled  airframe  dynamics  modifies 
the  effective  pitch  attitude  numerator,  1^^,  as  shown  in  Eq.  2. 


<p-^6a 


9 9 9 

Ngg  + KqjI^aBe 


= -2.9(-.5)(.86)[.28,  .i8][.27,  1.U]  (2) 

+ K,p|-i.5(-.OU6)(.29)[.5,  .6]f 


That  is,  the  effective  9 6^  numerator  is  the  sum  of  the  basic  open-loop 

0 ^0 

pitch  numerator,  %g,  and  a coupling  numerator,  N5j^6g,  t rm  "weighted"  by 
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the  roll  loop  gain,  K<p.  Ihe  roots  of  each  term  are  shown  in  Eq.  2 evaluated 
at  a = l8.8  deg  and  p = 6 deg.  If  Kq,  is  small  (i.e.,  rp  -*-Sa  essentially 
open  loop),  then  the  roots  of  are  essentially  the  same  as  those  of  N§g. 
However,  if  Kq,  is  very  large  then  the  coupling  term  dominates  and,  within 
the  frequency  range  of  interest,  the  roots  of  N§g  will  be  those  of  N^aig. 

For  intermediate  values  of  Kqj  the  root  l.''''ations  may  be  defined  in  terms  of 
the  root  locus  1 + (Kq^j^fg/l^g),  which  is  sketched  in  Fig.  8.  Note  that 


9 * 

Figure  8.  Ngg  Root  Migration  with  Increasing  Kq> 

both  the  numerator  and  denominator  have  roots  in  the  right  half  plane. 
Therefore,  1 (the  root  emanating  from  l/Te^)  is  driven  toward  and  even- 
tually into  l/T0r^,p  the  coupling  numerator  term  as  Kq)  increases.  Thus, 
closure  of  the  cp  6^  ]oop  decreases  the  potential  divergence  rate  caused 
by  the  9 ♦ 6g  closure  but  is  not  able  to  completely  stabilize  the  motion 
because  a zero  of  Ngg  always  remains  in  the  right  half  plane.  In  fact,  the 
first-order  divergence  will  lead  to  eventual  saturation  of  aileron  control 
(if  used)  as  the  divergence  progresses. 
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xt  is  important  to  recognize  that  in  an  analogous  manner  to  Fig.  7 
closure  of  the  0 6e  loop  will  modify  the  <p  -**•  6a  numerator  as  shown  in 
Eq.  5. 


= .UU[.15,  .i6][.55,  1.9]  (5) 

+ KeH.5(-.Oi^6)(.29)[.5,  .6]| 


Note  that  both  Eqs.  3 and  2 have  the  same  coupling  numerator.  Therefore, 

<p« 

tight  closure  of  the  0 6e  loop  will  result  in  the  effective  numerator  N5a 
having  a right  half  plane  zero.  Tight  closure  of  both  q)  6a  and  0 ♦ 6e 
will  result  in  identical  effective  numerators  having  right -half -plane  zeros 
in  both  loops  and,  since  the  closed-loop  poles  are  driven  into  the  zeros,  the 
same  divergence  in  each  loop. 

The  sensitivity  of  the  right-half -plane  zeros  l/Te^  suid  l/Teqp^  to  change 
in  trim  sideslip,  angle  of  attack,  or  bank  angle  is  shown  in  Table  3.  It  can 
be  observed  that  increasing  angle  of  attack  or  sideslip  (or  both)  rapidly  moves 
both  zeros  further  into  the  ri^t  half  plane,  and  divergence  time  constants  of 
two  seconds  or  less  would  be  expected. 


TABLE  3 

EXAMPLE  VARIATION  OF  RIGHT-HALF-PLANE  ZEROS 
DUE  TO  PITCH-ROLL  COUPLING 


Oq 

Po 

i/Tej 

1 /l0cpi 

18.8 

6 

5.11+ 

-0.3 

-0.05 

i8.8 

6 

0 

-0.U5 

-0.15 

20.9 

6 

5.  i 

-0.1+42 

-0.41+5 

19.5 

13 

0 

-0.72 

-0.31 

2k 


If  the  pilot  were  to  use  rudder  instead  of  aileron  to  oppose  the  lateral 
divergence,  the  results  would  be  similar  to  those  above,  because  the  coupling 
numerators  and  sN§g|j.  also  have  a ri^t-half- plane  root: 


»5e?r 

= -5(-.0l6)(.53)(2.05)(-2.k7) 

('■*) 

snlJr 

= 5.1(  .0lk)[.99,  .5.^][.25,  1.25] 

(5) 

Wius,  use  of  rudder  along  with  elevat  >r  will  not  stop  the  divergence  and  will 
result  in  saturation  of  the  rudder  control  as  well. 

Q.  Model  Validation 


It  was  previously  indicated  (Fig.  5)  that  a 5 model  of  the  pitch 
transfer  function  exhibit  the  same  key  pole-zero  values  as  does  the  6 DOF 
model.  Figure  9 snows  time  traces  of  a departure  obtained  witr  a nonlinear 

5 DOF  analog  simulation  of  the  A-J.  As  in  Fig.  5,  the  X equation,  or  speed 

degree  of  freedom,  has  been  eliminated.  The  same  trim  conditions  prevail 
as  for  the  6 DOF  time  traces  of  Fig.  6:  Oq  = l8.8  deg,  8o  = 6 deg,  level 

flight,  pitch  attitude  loop  closed  with  unity  dc  gain,  and  a 0.01  i*ad  9c 
introduced  to  disturb  the  trim.  The  dashed  lines  on  the  a,  p,  and  p traces 
indicate  the  response  obtained  from  the  nonlinear  6 DOF  digital  model.  Hie 

6 DOF  and  5 20F  time  traces  are  idertical  for  yaw  rate.  The  dotted  curves 
in  Fig.  9 stre  time  responses  for  a,  p,  P,  and  r obtained  from  an  actual  nose 
slice  departure*  in  an  A-7  aircraft  (Ref.  17)  but  with  the  flight  trace  shifted 
in  time  to  achieve  an  effect... ve  overlay  of  p traces.  The  close  agreement 
between  the  three  sets  of  time  traces  is  considered  to  be  a validation  of 

both  the  five-  and  six-degree-of-freedom  models. 

In  Fig.  10  the  m ♦ &a  loop  has  been  closed  along  with  9 &e*  This 

results  in  somewhat  closer  agreement  between  flight  and  the  5 ^F  analog  model 
a and  p traces.  However,  the  sideclip  is  not  in  quite  as  close  agreement  as 
it  is  in  Fig.  9.  It  was  subsequently  learned  (Ref.  18)  that  the  flight  traces 
reflect  use  of  rudder,  not  aileron,  to  oppose  the  norc  slice.  The  rudder  was 

*W  = 22,699  lb  Xcg  = 29. b i MAC 

h = 1k,l+00  ft  = 75  deg 
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Figure  10.  Departure  Time  Traces  with  Both  Pitch  and  Roll  Loops  Closed 


iiieffective  in  preventing  the  departure,  and  showed  a first-order  divergence 
to  full  deflection  similar  to  that  shorn  for  the  aileron  trace  in  Fig.  10. 

This  is  consistent  with  the  previous  coupling  numerator  analysis.  Also, 
coo^iarlng  the  time  traces  of  Fig.  9 and  10,  closure  of  the  bank-angle -to- 
alleron  loop  reduces  the  departuire  rate,  but  does  not  prevent  departure} 
verifying  the  analytic  prediction. 

Several  euMltional  points  emerge  flrom  an  examination  of  Fig.  10  as 
follows.  The  roll  rate  remains  relatively  small  until  approximately  12.5  sec 
beyond  ^diich  it  builds  very  rapidly  for  both  the  five-degree-of-fJreedom  model 
and  for  the  actual  aircraft.  !Qie  yaw  rate  traces  match  pretty  well,  but  the 
coopted  sideslip  is  approximately  twice  that  for  the  actual  aircraft.  Angle 
of  attack  for  both  fli^t  and  coo^nited  traces  match  very  well.  Thus,  based 
on  either  coopted  or  actual  time  histories,  the  aircraft  motion  initially  is 
a yawing,  sideslipping  response  which  rapidly  builds  to  large  values  before 
an  appreciable  roll  rate  develops  (as  described  in  Ref.  17). 

Relative  to  stall  and  post-stedl  recovery.  Refs.  I7  and  18  indicate  that 
such  is  best  accon^jlished  by  neutralizing  the  cockpit  controls  or,  preferably, 
"letting  go"  of  the  stick.  This  opens  the  loop  and  returns  the  aircraft  roots 
to  the  stable  open-loop  pole  positions,  as  shown  in  Fig.  U}  the  aircraft  then 
should  recover  by  itself. 

On  the  basis  of  the  foregoing,  it  is  concluded  that  the  right -half -plane 
0 

zero  of  the  Ngg  numerator  is  the  key  parameter  in  the  A-7  nose  slice  departure. 
Tite  relatively  close  cgreement  among  analysis,  simulation  traces,  actual 
flight  traces,  and  recovery  procedures  recoonended  by  LTV  is  considered  to 
validate  our  model. 

3 . Source  of  /a,  Na 

Preceding  analysis  has  shown  that  the  zero  I/T95  is  strongly  influenced 
by  the  stability  derivatives  and  Hi.  These  are  generated  in  the  9x9  matrix 
of  Fig.  2 as  the  partial  derivatives  of  p with  respect  to  a and  of  r with  respect 
to  a,  as  shown  in  Table  4.  These  partieds  are  expanded  in  both  literal  and 
numerical  terms  in  Table  4.  In  both  cases  the  first  term  in  the  sum  of  peirtial 
derivatives  dominates:  jCa  from  the  x>artial  of  roll  moment  with  respect  to  a 
but  at  a trim  6e  and  Pj  from  the  partial  of  the  yaw  moment  with  respect  to  a 

26 


but  at  a trim  6g  and  p.  Plots  of  these  coefficients  are  presented  in  Figs,  11 
and  12.  Figure  11  shows  Ct  is  positive  over  the  region  15-22  deg  a.  At 
approximately  22-23  deg  n.  the  slope  is  zero,  and  for  higher  angles  of  attack 
the  slope  reverses.  The  slope  also  increases  as  trim  sideslip  increases. 

Figure  12  shows  a decided  break  in  Cuq,  at  approximately  15  deg  a and  large 
negative  slope  to  approximately  5^  deg  where  the  slope  becomes  zero.  For 
larger  a this  coefficient  becomes  positive.  Thus,  in  the  region  between 
13  and  2h  deg  a,  there  is  a large  yawing  moment  with  change  in  ar.gle  of 
attack. 

Referring  back  to  the  time  traces  of  Figs.  9 or  10  and  the  matrix  of 
Fig.  2,  a strong  correlation  between  the  aircraft  motions  and  the  slopes  C/q, 
and  Cna  can  be  detected.  For  example.  Fig.  2 indicates  that  at  roughly  19  deg 
a and  6 deg  e,j(i  and  are  of  opposite  sign  and  of  nearly  the  same  magnitude. 
Thus,  with  p and  a increasing  positively  we  would  expect  the  resulting  rolling 
moments  contributions  to  nearly  cancel.  This  is  reflected  in  both  Figs.  9 and 
10  as  a relatively  low  roll  rate  up  to  an  angle  of  attack  of  approximately 
21-22  deg.  Above  this,  changes  sign,  X<i,  andjCp  then  augment  one  another 
and,  on  the  time  traces,  roll  acceleration  suddenly  becomes  very  large.  The 
large  negative  value  of  Mi  results  in  a large  negative  f and  positive  p as  a 
is  increased.  Again,  this  is  reflected  in  the  time  traces  of  Figs.  9 and  10. 
Thus,  the  departure  characteristics  of  the  A-7  appear  closely  related  to  the 
Cho,  and  C£^  characteristics  between  15  and  2k  deg  a in  Figs.  11  and  12.  Similar 
conclusions  are  reported  in  Ref.  19  regarding  the  effect  of  Cna  and  at 
P ^ 0 (e.g..  Fig.  15)  influencing  the  spin  entry  and  recovery  characteristics 
of  one  version  of  the  F-100. 

Results  of  wind  tunnel  tests  conducted  by  LTV  are  reported  in  Ref.  I8.  Of 
particular  interest  are  tuft  studies  showing  airflow  over  the  wing,  aft  fuse- 
lage, and  vertical  tail  at  16-19  deg  a and  O-16  deg  p.  Figure  lL  (from  Ref.  I8) 
depicts  tuft  activity  over  the  wings  at  17  deg  a and  various  sideslip  angles. 

At  17  deg  a and  zero  p,  an  area  of  steady  flow  exists  over  a portion  of  the 
center  section  of  each  wing.  As  p increases  the  area  of  sepeuration  moves 
rapidly  inboard  on  the  leading  and  outboard  on  the  trailing  wing  center 
sections.  The  larger  area  of  separated  flow  on  the  leading  wing  does  not  mean 
that  surface  will  stall  first,  since  tufts  are  not  necessarily  reliable 
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Figure  11.  Cj^ia,  p,  &g)  vs.  a at  6g  = -15' 


HCGATIVE 


-lOOr'  Static  Moment  Coefficients  Used  in  IBM  Calculations  (from  Ref.  19 


Figure  llj-.  A-7  Wing  Tuft  Activity;  Basic  Airplane  (from  Ref.  lB) 
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indicators  of  the  level  of  lifting  performance.  However,  there  is  strong 
evidence  that  separated  flow  develops  on  the  wing  center  section  coincident 
with  the  a's  and  P's  which  show  a marked  reduction  in  airplane  directional 
stability.  Figure  15  (also  from  Ref.  18),  for  19  deg  a and  6 deg  p,  shows 
a disturbance  of  the  tufts  near  the  tip  and  on  the  downwind  or  lee-side  of 
the  fin.  nie  orientation  of  the  tufts  on  the  fuselage  (not  shown  in  Fig.  15) 
and  on  the  vertical  tail  indicate  the  disturbemce  is  due  to  a vortex  system 
which  originates  at  the  intersection  of  the  wing  and  fuselage.  It  would 
appear  that  if  sideslip  were  increased  to  greater  than  6 deg  this  vortex 
would  clear  the  vertical  tail  and  stability  would  again  increase.  However, 
the  LTV  studies  showed  the  pattern  of  wing  center  section  flow  separation 
with  a and  P is  such  that  the  separation  on  the  upwind  panel  and  the  attached 
flow  on  the  downwind  panel  induces  a vortex  system  off  the  wing  surface  >diich 
has  the  same  sense  of  rotation  as  the  wing  fuselsige  vortex  just  noted.  Thus, 
the  two  systems  reinforce  each  other  and  create  a strong  vortex  which  results 
in  the  local  flow  redirection  at  the  vertical  tail  leading  edge.  This  flow 
orientation  would  produce  destabilizing  forces  and  account  for  the  directional 
stability  decrease  with  increasing  a and/or  p.  The  LTV  tuft  studies  indicated 
that  both  the  leading  and  trailing  wings  stalled  at  between  2h  and  25  deg  a. 

This  is  quite  consistent  with  the  zero  slope  points  on  Fig.  li. 

Ihus,  Ref.  18  relates  the  decrease  in  directional  stability  to  the  vortex 
activi+'y  from  the  fuselage  and  wing  center  panel  impinging  on  the  downwind 
side  of  the  vertical  fin.  It  appears  that  this  activity  can  equally  be  related 
to  the  inflections  in  the  coefficient  plots  of  Figs . 1 1 and  1 2,  to  the  sta- 
bility derivatives /a  and  Na  and,  eventually,  to  the  nose  slice  parameter,  T/Tq-^. 

k,  Concluilons 

The  foregoing  analysis  and  simixlation  has  indicated  that  the  nose  slice 
departure  of  the  A-7  aircraft  can  occur  at  a’s  considerably  below  the  normal 
stall.  The  nose  slice  is  caused  by  a chain  of  events  which  evolve  from 
directional  mistrim  or  failure  to  suppress  sideslip  in  maneuvers.  This 
sideslip  results  in  shed  vortices  which  give  rise  to  aerodynamic  moments, 

Na  andXa.  A nonlinear  kinematic  coupling  also  occurs  between  the  sideslip 
and  yawing  or  rolling  of  the  aircraft  to  give  rise  to  pseudo-derivatives,  Zp 
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and  Zy  (see  Fig.  2).  These  combined  coupling  effects  result  in  a right -half- 
plane zero  in  the  pitch-attitude-x.o-elevator  numerator.  Normal  pitch  attitude 
control  activity  then  drives  the  aircraft  unstable.  The  instability  starts 
as  a first-order  divergence  in  lateral,  directional,  and  longitudinal  modes. 
Pilot  atten^jts  to  oppose  the  motion  rapidly  leswi  to  t>aturation  of  control  in 
all  axes.  The  principal  meeuis  of  recovery  is  to  open  the  pitch  attitude  loop, 
that  is,  to  let  go  of  the  stick. 

A quasi-linearized  five-degree-of-freedom  model  containing  the  coupling 
terms  Xa,  Na,  Zp,  and  Zy  only,  and  with  appropriate  derivatives  changed  as  a 
fu;.''%ion  of  angle  of  attack  and  sideslip  provides  very  good  agreement  with 
actual  flight  traces.  This  model  can  be  utilized  for  the  piloted  simulations 
discussed  later. 

B.  wm  ROCK 

Another  common  departure,  or  pre -departure,  characteristic  is  wing  rock. 

This  generally  is  described  as  a fairly  large-angle  roll  oscillation  but  there 
is  considerable  disagreement  as  to  whether  it  is  an  open-loop  or  closed- loop 
('i.e.,  pilot-aggravated)  phenomenon.  Wing  rock  is  exhibited  by  several  current 
high-performance  aircraft j however,  the  F-I4  has  received  the  most  notoriety  in 
this  respect.  Wing  rock  reportedly  is  encountered  in  the  F-4  starting  at  about 
2t  deg  a and,  if  a is  further  increased,  this  is  followed  by  nose  slice  at  about 
2^  leg.  The  underlying  open-loor  cause  has  been  attributed  to  nonlinear  Cnp 
'Refs.  J and  20)  and 'or  the  partitioning  of  the  total  roll  damping  coefficients 

Cfr  - cos  a 
C£p  — C£p  sin  a 

into  the  basic  p,  r,  and  p contributions  (Ref.  21).  The  underlying  closed-loop 
cause  is  generally  attributed  to  "rlassical"  ‘-Aid  effects  fe.g..  Ref.  6). 

Several  F-4  siraulations  of  varying  degrees  of  complexity  and  fidelity  have 
been  made  to  exhibit  "wing  rock"j  however,  the  principal  means  of  achieving 
it  have  also  varied  considerably.  For  example.  Ref.  22  indicates  wing  rock 
could  only  be  obtained  when  ail  of  the  following  were  present  simultaneously. 
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1 . With  sane  form  of  control  (pilot  or  autc^ilot ) to 
prevent  the  aircraft  from  departing  before  wing 
rock  could  develop. 

2.  Ety  assuming  an  assmnetrical  moment  frcm  the 

stabilator  (i.e.,  and  N&g). 

3.  With  cp  6a  loop  closed. 

The  asynmetrical  joment  from  the  stabilator  was  required  to  start  a lateral 
motion.  It  should  be  noted  that  this  study  did  uct  Include  an  coupling 

derivative.  Mp  will  be  shown  lat<  to  have  a significant  magnitude  for 
the  Another  approach  to  achieving  closed-loop  wing  rock  has  been  to 

significantly  reduce  the  damping  derivative  Lp. 

In  an  effort  to  gain  better  insight  and  physics^,  under stauiding  of  this 
phenomenon,  the  tack  taken  here  has  been  to  review  the  potential  causal  factors 
noted  previously,  to  seek  new  parameters  throu^  the  same  analytic  techniques 
applied  to  nose  slice,  and  finally,  to  incorporate  the  most  promising  of  these 
causal  factors  into  th?  A-7  dynamic  model  to  see  if  wing  rock  can  be  induced 
in  this  aircraft  which  has  never  been  reported  to  have  a wing  rock. 

1 . Dynamic  Model 

Nonlinear  aerodynamic  data  for  the  F-1  were  supplied  in  the  form  of  digital 
lookup  tables  as  f(a,  g)  in  listings  and  card  decks  by  Mr.  J.  R.  Chambers  of 
NASA  Langley  Research  Center.  The  trim  configuration  selected  was  based  on  a 
flight  condition  and  loading  which  produced  wing  rock  during  the  flight  tests 
reported  in  Ref.  25." 


ii  = 14,000  ft 

Ix  = 26,000 

W = 57,200  lb 

ly  — 1 66,000 

c.g.  = 29.3  C 

Iz  = 185,500 

Flaps  and  gear  up 

^xz  = 5,000 

Open-loop  time  traces  from  the  nonlinear  digital  simulation  triumed  for  level 
flight  at  a = 20  deg  and  g = 1 deg  are  shown  in  Fig.  16.  A negative  ruddar 
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pulse  «a»  used  to  excite  the  unaugsiented  airframe.  Althou^  the  traces 
eidhibit  an  oscillation  that  is  obviously  self-sustaining  (see  p trace)  it 
bears  little  resemblance  to  the  large  amplitude  roll  usually  seen  in  flight 
traces  and  therefore  cannot  be  considered  wing  rock  per  se. 

Stability  derivatives  calculated  for  three  different  combinations 
of  trim  a and  3 are  presented  in  Table  5.  One  case  for  the  A-7  is  also 
3hown  for  con^ierlson.  Derivatives  reflecting  significant  differences  between 
the  two  aircraft  are/p,  Xr,  1^,  Hjp,  and  Of  these  the  most  startling  are 
jCp,  idiich  has  a change  in  sign  and  an  order- of -magnitude  difference,  and  Mp, 
which  has  no  counterpart  in  the  A-7  data.  Dote  that  and  Xa  are  comparable 
for  similar  trim  conditions. 

A plot  of  C|n  versus  p is  shown  in  Pig.  17^  where  it  can  be  seen  than  Cm^ 
is  negligible  for  p < deg  and  a < 15  deg  bi’t  becomes  significant  for 
a = 17.5  deg  and  above.  This  is  precisely  the  region  where,  wing  rock  is 
encountered.  A check  of  the  1^.  contribution  to  vehicle  dynamics  was  made  by 
calculating  the  transfer  functions  with  and  without  Ihe  results  are  shown 
in  Fig.  18  for  the  = 20  deg,  = 1 deg  case.  It  appears  that  for  this 
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Figure  17.  Cm  Versus  p from  F-1  Lookup  Tables 


TABLE  S.  BODY  AXIS  DIM0ISIONAL  DERIVATIVES  (F-4  AND  A-7) 
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Figure  i8.  Effect  of  Mp  on  F-U  Poles  and  Zeros 

(oo  = 20C,  Po  = lO) 
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case  the  principally  affects  the  lateral  and  longitudinal  short-period 
inodes  with  a resultant  increase  in  ^gp  and  decrease  in  Further,  it 
appeaurs  that  with  =-1.63,  closure  of  any  or  all  of  the  three  loops 
shown  in  Fig.  l8  could  easily  result  in  an  unstable  dutch  roll  mode.  These 
loops  could  be  closed  by  the  pilot  (0  and  ip)  ani  the  yaw  rate  damper. 

The  effect  of  was  verified  with  a 5 DOF  analog  simulation  which  pro- 

duced a stable  dutch  roll  for  M(-  = 0 and  divergent  dutch  roll  of  =—1.6. 

[it  should  be  noted  in  passing  that  a similar  shift  in  short-period  modes 
was  obtained  in  the  Ref.  l6  analysis  of  nonsynmetric  flight  effects.  However, 
this  was  attributed  to  the  coupling  derivative  /i. ) 

Figure  17  indicates  that  Cm^  changes  sign  depending  upon  the  sign  of  p 
such  that  is  always  negative.  For  positive,  linear  analysis  indicates 
the  migration  of  the  longitudinal  short-period  and  dutch  roll  modes  to  be 
reversed  from  that  shown  in  Fig.  18.  That  is,  the  longitudinaJ.  short  period 
moves  toward  the  right  half  plane  while  the  dutch  roll  moves  to  the  left. 
However,  as  shown  below,  other  factors  in  the  nonlinear  6 DOF  model  apparently 
bias  the  aircraft  motion  towards  an  unstable  dutch  roll. 

At  the  "equest  of  the  USAF,  McDonnell-Douglas  Aircraft  Co.  provided  a 
collection  of  wing  rock  traces  obtained  during  flight  test  of  the  slatted- 
wing  aircraft.  TypiceJ.  of  these  are  the  1 g stall  traces  shown  in  Fig.  19 
(reproduced  from  Ref.  2U)  which  reflect  trim  conditions  very  similar  to  those 
selected  for  the  preceding  analysis,  i.e.: 

h - 17,500  ft 
W = 57,^»10  lb 

c.g.  = 29.5^  M^C 

M = 0.56 
Clean  aircraft 


The  traces  indicate  that:  ' 

a.  Wing  rock  is  not  a controls-fixed  oscillation  but 
appears  to  be  caused  by  or  aggravated  by  the  pilot 
(i.e.,  a closed- loop  phenomenon). 


b.  It  is  not  uncoonon  to  have  a bias  sideslip  of  one 
or  acre  degrees  during  wing  rock. 

c.  Relatively  large  bank>angle  and  sideslip  excursions 
but  SBiall  angle-of-attack  oscillations  are  exhibited. 

It  was  determined  from  the  traces  that  the  average  gain  between  roll  angle 
and  lateral  stick  excursions  reflect  a pilot  gain  Kp^  = 0.2  deg/deg.  Ihere- 
fore,  a closed- lo<^  simulation  was  run  utilizing  the  nonlinear  6 1X)F  airframe 
and  a sln5)le  pilot  model,  Yp  = 0.2  e~*^®.  The  resulting  traces  are  shown  in 
Fig.  20  for  an  initial  = -t-1  deg  and  Pig.  21  for  an  initial  fo  = ~1  deg. 

These  attempts  at  modeling  the  real  aircraft  traces  show,  in  Fig.  20,  that 
roll  oscillations  for  ^o  = +1  deg  grow  much  greater  than  the  longitudinal 
oscillations.  The  roll  divergence  has  approximately  a 10  sec  time  to  double 
aii5)litude  €Uid  hence,  over  a time  i>eriod  of  50  to  UO  seconds,  coiild  develop 
into  a roll  rate  of  about  50  deg/sec  amplitude  which  would  approximate  that 
shown  in  'Me  fli^t  trace  of  Fig.  19.  Figure  21  (Po  = deg)  does  not  show 
the  same  oscillating  divergence  tendencies  but  it  should  be  noted  that  the 
sideslip  is  steadily  building  toward  a i>ositive  value  which  would  then  diverge 
as  in  Fig.  20.  This  is  particularly  interesting  since  the  majority  of  the 
flight  traces  in  Ref.  24  indicate  the  same  bias  toward  AHL  (positive)  p.  Thus, 
it  appears  that  positive  p,  negative  and  pilot  closure  (intentional  or 
Inadvertent)  of  roll  (attitude,  rate,  or  acceleration)  to  aileron  may  be  a 
significaiit  contributor  to  wing  rock. 

Another  investigation  of  possible  sources  of  the  F-4  wing  rock  included 
separation  of  the  total  roll  damping  coefficients 

Cjjj.  Cfp  cos  a 
Cfp  + sin  a 

into  the  basic  p,  r,  euid  p contributions.  The  separation  was  based  upon  the 
coefficient  data  of  Ref.  :^1  and  resulted  in  the  following  sets  of  derivatives 
for  the  Oq  = 25  deg,  Pq  = 0 deg  case: 
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Note  this  is  for  symmetric  (Pq  =0)  flight.  The  effects  were  assessed  via 
litersJ.  e3q>emsion  of  the  three-degree-of-freedom  lateral  matrix  and  an  analog 
simulation.  Both  showed  little  difference  in  the  resulting  lateral-directional 
dynamic  characteristics  and  verified  the  acceptability  of  lun5)ing  all  roll 
damping  into  an  effective  Xj  and 

The  possible  influence  of  nonlinesur  Cn^  (Refs.  5 and  20)  was  also  inves- 
tigated using  the  five-degree-of- freedom  analog  simulation.  In  this  instance 
the  nonlinearity  in  Np  was  modeled  after  the  Cnp  slopes  indicated  in  Refs.  3 
and  20,  i.e.,  reverseJ-s  at  5 to  10  deg  p.  The  results  were  similar  to 
those  of  a bang-bang  control  system  operating  with  a large  dead  zonej  the 
aircraft  drifted  through  the  unstable  p range  and  then  developed  a limit  cycle 
type  oscillation  about  the  stable  Cm  vs.  p intercept  point  (e.g.,  10  deg  p). 

This  is  not  compatible  with  fli^t  traces.  Similar  nonlineeir  representations 
of  CXp  and  Cnp  were  atten^jted  with  no  in5>rovement  in  results. 

Since  the  most  promising  causal  factors  appeared  to  be  the  basic  longitudinal- 
lateral  coupling  of  the  F-t  in  nonsymmetric  flight,  an  attempt  was  made  to  induce 
wing  rock  in  the  A-7  by  modifying  key  stability  derivatives.  The  principal 
difference  between  the  derivatives  of  the  F-l<  and  A-7  at  comparable  high  a has 
been  shown  previously  to  be  in  X^,  N^,  Xp>  Np,  and  Mp.  Therefore,  these  A-7 
derivatives  were  rescaled  on  the  basis  of  the  differences  in  vehicle  mass, 
inertia,  wing  area,  etc.,  and  new  derivatives  obtained  for  a "pseudo  F-4." 

These  are  summarized  as  follows: 
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A-7 

PSEUDO 

20 

18.8 

18.8 

1 

6 

6 

1.268 

0.3323 

1.21 

0.1906 

0.16 

K 

-0.505 

-i.»v25 

-0.67 

H 

0.52 

0 

0.424 

-1.63 

0 

-1.54 

These  new  values  were  incorporated  in  the  nonlinear  five-degree-of-freedom 
nnalng  simulation  of  the  A-7  and  evaluated  as  a potential  configuration  for 
the  piloted  simulation  discussed  in  Section  IV.  The  change  in  derivatives 
did  result  in  an  oscillatory  behavior  quite  similar  to  that  shown  in  Fig.  l6 
but  having  a lower  frequency  and  larger  junplitude  about  all  three  axes . When 
assessed  by  the  evaluation  pilots  in  a target  tracking  task  the  behavior  was 
not  considered  representative  of  wing  rock  (open  or  closed  loop).  Rather, 
the  nose  of  the  aircraft  tended  to  describe  a circular  motion  about  the  target. 
This  simulation  was  not  pursued  further. 

2.  Summary  and  Conclusions 

Various  open-  and  closed-loop  parameters  and  aerodynamic  coefficient 
anomalies  which  have  been  suggested  in  the  technical  literature  as  potential 
causal  and/or  contributing  factors  ;o  the  F-lt  wing  rock  have  been  investi- 
gated. The  only  situation  found  t'  produce  open-loop  oscillatory  character- 
istics in  a nonlinear,  six-degree- of- freedom  digital  model  is  nonsymmetric 
(p  / 0)  flight  at  high  c.  The  principal  factor  appears  to  oe  the  coupling 
derivative  which  has  a destabi..  Izing  influence  on  the  short-period  ^ lateral 
and/or  longitudinal)  modes.  The  resulting  oscillatory  motion  is  apparent  as 
angular  rates  and  displacements  about  all  three  aircraft  axes  of  approximately 
equal  magnitudes  and  therefore  does  not  specifically  resemble  wing  rocK. 

Study  of  all  available  flight  traces  which  present  aircraft  motion  and 
control  input  (or  surface)  provide  strong  evidence  that  the  F-t  phenomena  is 
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closed  loop.  Significant  laterad  stick  activity  is  generally  observed  to 
accon5)any  the  predominately  rolling/sideslipping  wing  rock  motion.  Incor- 
poration of  a simple  roll  control  pilot  model  in  the  above  nonlinear,  six- 
degree-of-freedom  digital  simulation  and  adjustment  of  the  feedback  gain  to 
match  that  observed  in  flight  traces  resulted  in  an  oscillatory  roll  diver- 
gence which  could  grow  into  wing  rock. 

An  aoten^t  to  change  the  departure  characteristic  of  the  A-7  aircraft 
from  nose  slice  to  wing  rock  proved  partially  successful.  Scaling  of  five 
key  stability  derivatives  to  resemble  those  of  the  F-4  did  result  in  an 
oscillatory  characteristic;  however,  it  did  not  have  the  predominate  rc.ll 
characteristic  of  wing  rock. 


SECTION  in 


SIMUIAnON  SESGRlPnai! 


This  section  describes  in  general  terms  the  five- degrees- of- freedon 
piloted  simulation.  Details  are  contained  in  Volucie  II,  Appendix  IV. 

The  simulation  involved  two  basic  tasks: 

• Assess  flying  quality  parameter  influences  on 

— Departure  susceptibility/resistance  (warning) 

- Departure  and  recovery  characteristics 

Pilot  flying  technique  in  the  stall/departure 
flight  regime 

• Measure  pilot  dynamic  control  peirameters  in  the 
staJl/departure  fli^t  regime 

To  provide  meaningful  data,  the  vehicle  dynamics  had  to  be  a realistic 
representation  of  hi^  angle- of- attack  flying  characteristics.  However, 
to  maintain  control  over,  and  vary  only,  specific  key  parameters  as  well 
as  to  avoid  unnecessary  ccmq)lication  which  might  confound  the  experiment, 
it  was  desirable  to  employ  the  simplest  possible  simulation  mechanization. 
As  indicated  in  the  previous  section,  a nonlinear  five-degrees-of-freedom 
analog  simulation  was  deemed  to  satisfy  both  sets  of  requirements.  The 
adequacy  of  this  simulation  was  validated  by  two  means.  First,  a dynamic 
check  was  made  aga-'nst  a nonlinear  six-degrees-of- freedon  digital  model. 
Secondly,  the  simulation  was  evaluated  by  two  U.  S.  Air  Force  experimental 
flight  test  pilots  who  are  experienced  in  the  A-7  aircraft.  One  pilot 
flew  the  A-7  between  sessions  in  the  simulator  to  specifically  assess  and 
conpare  the  simulated  and  actual  vehicle  characteristics.  Following  a 
few  minor  adjustments  both  pilots  pronounced  the  simulation  to  be  a good 
replica  of  the  A-7  characteristics.  Once  vsJ.idated,  the  appropriate 
flying  quality  parameters  were  then  varied  and  a series  of  experiments 
performed  which  will  be  covered  in  later  sections. 

A.  caaSERAL  LAYOUT  AND  OPERATION 

The  functional  layout  of  the  simulation  facility  is  shown  in  Fig.  22. 

A fixed  base  filter  cockpit  was  employed  wh:  ch  contained  a itiU  complement 
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of  fli^t  instruBent  displays  Including  an  all-attitude  ball.  An  out-the- 
window  display  wan  provided  in  front  of  the  windshield  via  a large  CRT  screen. 
This  resembled,  somewhat,  the  basic  H.U.D.  of  the  A-7*  Control  system  feel 
was  adjusted  to  match  the  A-7  breakout  and  force  gradients  as  closely  as 
possible.  A pedal  shaker  was  activated  as  the  aircraft  neared  stall.  Air- 
frame dynamics  and  display  geometry  were  modeled  on  two  EAI  251R  con^niters 
with  associated  nonlinear  racks.  Ttie  equation  flow  diagram  for  the  aircraft 
and  display  systems  is  presented  in  Fig.  23*  The  left  side  of  this  figure 
reflects  mechanization  of  the  airframe  equations.  The  Icwer  ri^t  side 
reflects  mechanization  of  the  target  aircraft  and  the  differential  motion 
between  target  and  the  controlled  aircraft.  Lateral  and  normal  acceleration 
as  sensed  by  cockpit  instruments  were  also  calculated.  The  upper  ri^t 
portion  of  Fig.  23  reflects  the  con5>lications  in  mechanizing  the  derivatives 
which  varied  as  functicxis  of  initial,  pertiorbation  and  total  a and  p. 

B.  AIRCRAFT  MODEL 

The  five-degrees-of- freedom  airframe  model  was  obtained  by  removing  the 
X equation  from  the  matrix  of  Fig.  2 and  varying  six  key  derivatives  in  the 
lateral  equations  as  a function  of  a and  p.  These  equations  also  reflect 
the  analog  mechanization  limitations  which  required  assuming  tan  P = sin  p. 

The  adequacy  of  these  approximations  were  verified  by  the  comparison  of  the 
nonlinear  1;  DOF  analog  and  DOF  digital  simulations  as  already  reviewed  in 
Section  II  (e.g..  Figs.  and  10).  Since  the  aircraft  stalls  at  approxi- 

mately 29  deg  a and  it  was  desired  to  start  the  simulation  runs  somewhat  below 
stall,  straight-line  fits  were  made  to  the  derivatives  only  for  the  a range 
17  to  23  deg,  and  these  were  then  extrapolated  to  30  deg  a.  Similarly,  the 
derivative  variation  with  sideslip  was  fitted  over  the  P range  0 to  ±8  deg 
and  then  extrapolated  to  larger  p.  The  two  derivatives  which  varied  con- 
tinuously with  a and  p are  shown  in  Fig.  2h . varied  linearly  with  a until 
well  past  the  stall  and  linearly  with  p to  4 deg  where  it  became  constant 
for  further  increases  in  p.  X,^  varied  linearly  with  p for  all  a but  not 
with  a until  approximately  21  deg.  Above  21  deg  it  rapidly  decreased  with 
further  a increase  and  changed  sign  at  approximately  23  deg  a.  This  was 
n.odeled  as  a constant  slope  with  p to  21  deg  a and  then  a decreasing  slope 
v/lth  further  a increase. 
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Figure  22,  Simulation  Block  Diagram 


th  a and  P 


Three  derivatives  were  varied  linearly  as  a function  of  a only  (Fig.  25). 
These  are  Yy,  and  N^.  The  latter  changes  sign  at  approximately  19  deg  a. 
Hi  was  varied  linearly  as  a function  of  P only.  The  remaining  thirteen 
stability  derivatives  were  maintained  fixed  throu^out  the  simulations  runs 
althouj^  different  initial  values  of  some  were  en5)lQyed  to  achieve  desired 
changes  in  pertinent  flying  quality  parameters.  Numerical  values  sure  given 
in  Appendix  D of  Voltone  II. 

The  derivative  mechanization  described  above  provided  a relatively 
unconqplicated  airframe  model  yet  afforded  the  time  trace  matches  reflected 
in  Figs.  9 and  10. 

The  contributions  of  the  pitch  and  roll  rate  command  augnentation  sind 
yaw  stability  augnentation  systems  were  not  modeled.  In  the  actual  aircraft 
the  roll  rate  CAS  is  automatically  faded  out  In  this  a range  because  it 
produces  pro-departure  surface  deflections.  Therefore,  it  need  not  be 
simulated.  Analysis  of  the  other  two  fixed  gain  systems  indicated  their 
contributions  to  be  insignificant  at  the  low  dynamic  pressure  of  the  con- 
ditions being  simulated. 

C.  DEFABTURE/RECOVERY  DISPLAY  AND  TASK 

A task  appi'oximating  an  air-to-air  tracking  tail- chase  was  presented  on 
the  head-up  display.  A view  which  shows  the  head-down  and  head-up  display 
is  shown  in  Fig,  26.  A symbol  representative  of  own-aircraft  centerline  v-ts 
fixed  at  the  center  of  the  CRT  face.  The  target  and  horizon  then  moved  with 
respect  to  the  own- aircraft  centerline  symbol.  Tlie  head- up  field  of  view 
was  ±24  deg  in  azimuth  and  ±l8  deg  in  elevation.  The  subject  aircraft  could 
be  maneuvered  in  bank,  pitch  and  heeiding  relative  to  the  target 

in  Fig.  25).  Sustained  attitude  deviations  would  result  in  latereQ  or 
vertical  displacement  with  respect  to  the  target.  Only  range  to  target 
remained  fixed  at  2000  ft.  The  subject  aircraft  pitch  and  roll  attitude 
could  be  discerned  with  respect  to  the  horizon  (|J>,  0),  while  the  target 
aircraft  in  view  provided  a reference  for  judging  heading  or  turn  rate  from 
the  head-up  display.  Without  the  target  aircraft,  the  pilot  cculd  not 
detect  nose  slice  or  onset  of  departure  until  a rapid  roll,  discernible 
from  motion  of  the  horizon,  would  develop.  (A  similar  A-7  simulation  reported 
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Figure 


Figure  26.  Outside  and  Inside  Display 


in  Ref.  2^  had  only  a head>doim  disx>lay  and  the  haslc  departure  was 
described  as  a roll  departure.)  In  actual  fll^t  the  pilot  has  two  cues. 

One  is  the  lateral  motioi  of  the  nose  of  the  aircraft  with  respect  to 
cloud  texture  or  a target  aircraft.  The  second  is  a rather  significant 
side  force  as  sideslip  and  lateral  acceleration  develop.  To  partly  con^ 
pensate  for  the  lack  of  lateral  acceleration  cue,  a replica  of  the  ball 
of  the  turn  and  slip  indicator  was  also  displayed  on  the  CRT  screen. 

This  proved  helpful  to  the  pilot  in  detecting  that  something  was  happen- 
ing or  about  to  happen,  because  he  could  detect  motion  of  the 
peripheral  vision.  However,  once  the  ball  was  pegged  at  either  side  of 
its  displacanent  it  was  of  little  value  as  a cue  except  when  the  aircraft 
was  nearly  under  control  following  a departurej  the  pilot  could  then  actually 
look  at  the  ball  and  determine  the  rudder  input  required  to  achieve  final 
coordination  and  recovery.  During  recovery  with  the  target  or  horizon  line 
not  within  view,  the  pilot  could  revert  attention  inside  the  cockpit  to 
the  attitude  ball  which  displayed  4,  0,  T. 

Pilot  comnents  on  the  display  were  that  it  was  quite  similar  to  flying 
tail  chase  using  the  actual  A-7  head-up  display.  As  indicated  previously, 
the  target  was  required  to  be  in  view  for  the  pilot  to  obtain  a nose  slice 
cue  at  high  angles  of  attack.  A continuous  digital  readout  of  an^e-of- 
attack  and  sideslip  was  displayed  above  the  head-up  display  to  help  the 
e:q)erlmenter  and  pilot  aneuLyze  departures,  control  functions  and  recoveries. 
These  were  most  helpful  to  the  experimenter  while  observing  pilot  control 
activity  and  apparent  detection  of  departure.  These  digital  displays  were 
also  of  benefit  to  the  pilot  in  analyzing  the  departure  characteristics 
and  the  effects  of  his  control  inputs  by  placing  the  computer  in  hold  and 
then  observing  sdl  of  the  motion  quantities  displayed. 

D.  PILOT  PARAMETER  MEASUREMEIITS 

Measurement  of  pilot  dynamic  parameters  was  accon5)lished  through  use 
of  a Describing  Function  Analyzer  (DFA)  as  described  in  detail  in  Refs.  26 
and  27.  Briefly,  the  CFA  is  a device  for  making  on-line  dynamic  response 
measurements  of  manual  control  systems.  It  generates  a sum-of  sinusoids 
input  forcing  function  that  is  used  to  excite  the  dynamic  system  under 
consideration;  and  computes,  on-line,  the  finite  Fourier  transform  of  a 
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i(t)  = E A^,  cos^cj^f  t<^') 

i-1 


The  real  and  imaginary  components  of  the  Fourier  transform  are  integrated 
during  a run  and  available  immediately  at  the  end  of  the  measurement  period. 
The  finite  Fourier  transform  of  error  signal  at  input  frequencies  is 


E(>k)  = He[E(jo)k)]  + jIm[E(  jcDk)] 

where  ^ ^tp 

ReE(j'jjj,)  ==  J dt 

''R  o 

o 

Off-line  calciilationc  are  then  required  to  obtain  the  error- to- input  describ- 
ing function 


^^iC'  j-'k'' 


R' j '>)  ^ Re[E(j'ik)]  + jIm[E(ju)k)] 

T'jxk)  Ak 


and  the  desired  open-loop  (syr.tem  output- to- error)  describing  function 


Yol^J'%)  = TT-r- — T “ ■■ 
Yie( J^k) 


To  obtain  pilot  describing  functions  in  multiaxis  tracking  situations, 
past  laboratory  procedures  have  involved  simultaneous  introduction  of  separate, 
uncorrelated,  forcing  functions  in  each  axis.  This  permits  direct  determination 


of  describing  functions  for  each  axis  as  well  as  for  cross-talk  or  noise 
between  axes.  Unfortunately  this  does  not  provide  a realistic  forcing 
function  f<T  sinmlation  of  air-to-air  tracking  where  between-axis  target 
motion  is  correlated.  That  is,  there  sure  certain  basic  relationships 
between  tsurget  bank,  pitch  and  hesiding  for  any  useful  aircraft  maneuver. 

It  is  not  likely,  for  example,  that  a tsurget  vrill  bsmk  into  a turn  and 
then  pitch  away  from  the  tiim  (i.e.,  negative  g turn).  Such  apparent 
target  motions  sure  rejected  by  a pilot  as  sigJ.t  noise  due  to  buffet, 
pipper  wander,  etc.  Thus,  to  provide  a realistic  input  which  e^qperienced 
filter  pilots  would  track,  it  was  necessary  to  introduce  the  forcing  func- 
tion input  to  one  axis  at  a time  (Fig,  27)  and  depend  upon  stationary  run- 
to-run  performance.  The  as^litude  and  frequency  of  the  five  sinusoids 
enqployed  in  the  forcing  functions  sure  shown  in  Fig.  28. 

The  measurement  of  pilot  describing  functions  with  the  DFA  requires  a 
linear,  stationary  control  process.  Therefore  the  airframe  dynamics  were 
"frozen"  for  the  data  runs.  This  was  accot5>lished  by  fixing  sdl  airframe 
stability  derivatives  at  the  value  appropriate  for  the  given  trim  condition 
and  aircraft  configuration. 
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Figure  27.  Tracking  Task  Lcx>p  Structures 


Figure  28.  Forcing  Function  Spectrum 


SECTZOR  IV 


nxim  wjst  BUUMkisR  assessmesii 

This  section  presents  an  analysis  of  quantitative  euid  qualitative  results 
obtained  in  the  five  degree- of- freedcm  piloted  simulation  of  vai’ied  departure/ 
recovery  flying  qusility  parameters-  Attention  was  concentrated  on  nose  slice 
departures  and  the  effect  of  key  derivatives  on  departure  and  recovery  charac- 
teristics. Emphasis  is  placed  on  physical  understanding  of  the  d^arture  and 
recovery  dynamic  characteristics,  i-e.,  a correlation  of  c^en-  and  closed-loop 
modes,  pilot  commenteiry,  and  strip  chart  recordings- 

It  is  shown  that  previously  suggested  lateral-directional  stability  and 
departure  parameters  do  not  correlate  with  the  departure  characteristics  and 
pilot  ratings  obtained  here.  The  only  parameter  which  appears  consistent  is 
the  coupling  zero  of  the  0 -*-6e  numerator.  The  results  obtained  are  con- 
founded by  a significant  difference  in  piloting  technique  between  the  two 
test  subjects. 

A.  FILOIZIK}  TASK 

The  basic  task  was  to  track  a target  which  led  the  pilot  into  a stall 
and/or  departure.  The  initial  condition  placed  the  aircraft  wings-level 
in-trail  to  the  target.  A ramp  target  flight  path  input  of  1 deg/sec  was 
commanded  to  obtain  a consistent  rate  of  stall  approach.  The  pilot  then 
followed  the  target  aircraft  flight  path  until  departure  motion  was  detected 
and/or  he  felt  he  could  no  longer  follow  the  target  aircraft  and  still 
achieve  a successful  recovery.  The  task  thus  represents  a straight-ahead, 
paced,  stall  approach.  Since  pilots  normally  evaluate  stall  characteristics 
of  new  aircraft  configurations  by  performing  wings-level  1-g  stalls,  this 
type  of  display  and  tracking  task  assisted  the  pilots  in  correlating  the 
simulator  displayed  motions  with  previous  flight  experience;  and  in  assign- 
ing Cooper-Harper  ratings  to  the  various  configurations  and  maneuvers.  A 
wind-up  turn  target  following  maneuver  was  also  simulated.  This,  however, 
resulted  in  an  identical  tracking  task  to  the  straight-ahead  wings  level 
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t«fk  except  that  the  target  aod  <nm-aircraft  vere  banked  with  respect  to  the 
horizon.  ISie  pilots  could  detect  no  basic  difference  in  the  tasks  or  in  the 
motions  and  therefore  the  wind-up  turn  was  discontinued. 

The  departure  cue  was  onset  of  uncoomanded  nose  yaw  with  respect  to  the 
teurget  aircraft  or  the  rea3ir.ation  by  the  pilot  that  the  controls  vere 
ineffectual  in  retvuming  his  aircraft  to  the  desired  fll^t  path.  Two  con- 
trol strategies  were  enqployed  following  detection  of  departure.  The  first 
was  to  apply  sustained  resistive  control,  that  is,  attenpt  to  maintain 
existing  pitch  attitude  at  departure  detection  and  then  to  apply  lateral- 
directional  controls  as  necessary  to  resist  the  departure.  The  second  was 
to  relax  longitudinal  control  to  the  zero-force  stick  trim  position  or  to 
apply  forward  stick  force.  Five  separate  recovery  techniques  were  investi- 
gated with  this  relaxed  longitudinal  stick  control: 

A.  Aileron  and  rudder  neutral 

B.  Aileron  only  to  regain  control 

C.  Rudder  only  to  regain  control 

D.  Aileron  and  rudder  to  regain  control 

E.  Same  as  D except reduced  by  a factor  of  four 

°r 

In  addition  to  strip  chart  recordings  of  the  various  aircraft  motions, 
pilot  commentary  was  recorded  continuously.  This  conmentary  included  a 
description  of  the  departure  warning  and  its  adequacy,  control  technique 
employed,  recoverability  of  the  aircraft,  and  the  Cooper-Harper  rating  of 
the  recovery  task.  The  edited  pilot  ccxnmentary  is  presented  in  Volume  II. 

The  plan  was  to  start  each  run  with  the  aircraft  at  an  initial  condition 
of  l8.8  deg  a and  6 deg  p.  This  resulted  in  such  a rapid  depeurture  onset 
that  the  pilots  could  not  react  in  time  to  delay  the  departure  or  accomplish 
recovery.  The  sideslip  initial  condition  was  then  reduced  to  3 deg.  This 
proved  satisfactory  from  a control  standpoint;  however,  f.t  soon  became  appa- 
rent that  the  pilot  would  get  in  a "groove"  because  the  nose  slice  was  always 
in  the  same  direction.  Thus,  the  pilot  would  anticipate  the  nose  slice  motion 
and  apply  corrective  controls  in  an  essentially  precognitive  manner.  For  the 
final  data  runs  the  aircraft  initial  condition  was  set  to  18.8  deg  a and 
0 deg  p.  With  these  trim  conditions  the  aircraft  would  depart  in  either 
direction  and  this  depended,  in  part,  on  the  random  buildup  of  sideslip  due 
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to  control  activity  in  tracking  the  target.  This  initial  condition  also 
afforded  several  seconds  of  relatively  "precise"  tracking  and  allowed  the 
pilot  to  obtain  a feel  for  the  aircraft  as  it  approached  the  stall  and 
departure. 

The  analysis  to  follow  is  based  on  a representative  3 <isg  p since  this 
generally  reflects  the  initial  departure  situation  (i.e.,  warning  and 
detection) . 

E.  EIRANRBt  MAUAimS  AID  H8DLT8 

Nine  aircraft  dynamic  configurations  were  employed  in  the  simulation. 
Configuration  1 represented  the  nominal.  A-7  dynamics.  Key  stability  deriva- 
tives were  changed  (in  the  simplified  manner  described  in  Section  III-B)  to 
selectively  control  the  initial  condition  values  of  the  open-  and  closed- 
loop  handling  quality  parameters  l/T^,  oBR,  oq),  top/uB>  and  which 

then  varied  as  a function  of  a and/or  ^ as  the  aircraft  approached  stall.  The 
strategy  was  to  first  worsen  the  coupling  divergence  parameter  (1/T05)  and 
selectively  Improve  the  lateral  handling  parameters  as  indicated  for  Configu- 
rations 1 throu^i  3 in  Table  6.  Secondly,  the  divergence  parameter  was 
improved  and,  again,  the  individual  lateral  parameters  were  varied  as  indi- 
cated for  Configurations  4 through  9 i*i  Table  6.  Representative  values  of 
the  parameters  (and  derivatives)  at  = 18.8  deg  and  = 3 deg  for  all 
nine  configurations  are  tabulated  in  Table  7. 

1 . CoDflguxmtion  1 

Figure  29  shows  a typical  migration  of  airframe  open-loop  poles  and  zeros 
with  fixed  a and  varying  p.  The  arrows  indicate  movement  of  poles  or  zeros 
as  p is  increased.  The  inportant  points  are: 

• There  is  relatively  slight  movement  of  all  pOles 
for  0 < p < 6 deg;  dutch  roll  damping  decreases. 

0 

• The  two  N5  real  zeros  move  toward  the  ri^t  half 
plane  and  there  is  large  increase  in  with  large 
increase  in  mid- frequency  anplitude, 

• 

pensation  for  signi^'^ant  closed-loop  control. 

^)7 


The  mid-frequency  zero  oj(p  moves  with  so  that 
cu^oB  remains  relatively  constant;  the  roll  trans- 
fer function  is  K/s^-like  and  requires  pilot  com- 


mam  6.  CONFlCiUKATION  CHAKES 


CBARACTEBISnC  CB 


VoHiiial  A-7 


locrcaae  divergence  nte  of  aoee  eliee  aaaoelated 
with  KHP  zero  l/Te^ 

BHalnate  edverse  eileron  yaw  and  aubaoqnaat  aaeitatioit 
of  Otttdi  roll  node 

Betain  coupled  spiral/roll  aiibsldenee  or  lateral 
phugoid  node 


Increase  divergence  rate  of  nose  slice 
Elininate  lateral  phugoid 

Moderate  decrease  In  aileron  excitation  of  Oatdi 
roll  node 


Ellalnate  nose  slice  associated  vlth  BHP  zero  l/Te| 

EliBinate  adverse  ail'sron  yaw  and  excitation  of  Dutch 
roll  mode 

Setain  lateral  phugoid 


Eliminate  nose  slice  associated  with  RHP  zero  l/Ie| 

Retain  large  aileron  excitation  of  Dutch  roll  and  near 
roll-reversal  characteristic  for  aileron  control  of  bank 

Retain  lateral  phugoid 

Increase 


Eliminate  nose  slice  associated  with  l/Te> 


Retain  near  noainal  u>)/<U(j 
Eliminate  lateral  phugoid 
Decrease  Kt^dyn 


Same  as  6 abc  v ■? 

Decrease  to  determine  effect  on  roll  rate  obtained 
in  departure  anu/or  recovery 


Provide  RHP  zeros  in  both  longitudinal  and  laterad 
attitude  control  tasKs 

Provide  initial  RHP  poles  (lateral  phugoid) 

Both  of  above  acccmplisned  by  making  %<iyn  ” ® 

Provide  nuniua'  nice  slice  associated  with  l/Tflj 

Provide  "gool"  roll-subsidence  and  spired  cheiracteristlcs 


Degraic  Ltit?;;  to.:  damping 


L 


MATRIX  OF  IIANDLING  PARAMETER  VALUES 


• The  total  absence  of  mid- frequency  zeros  and  the  presence 
of  a ri^t  half  plane  zero  in  result  in  large  low- 
and  mid- frequency  phase  lag  and  make  this  control  tech- 
nique inq)ractical. 

• The  migration  of  the  hi^-frequency  zeros  of 
increases  separation  from  the  short  period  poles  and 
thereby  indicates  increasing  excitation  of  the .longi- 
tudinal mode  if  the  pilot  attempts  to  control  ♦ with 
aileron. 

'|f 

• The  movement  of  the  mid-frequency  zeros  of  Nsp  to 
stay  close  to  the  dutch  roll  poles  indicates  little 
effectiveness  of  the  rudder  in  damping  dutch  roll; 
however,  the  transfer  function  remains  K/s-like  and 
provides  relatively  good  low-frequency  control  of 
yaw  without  pilot  compensation. 

Figure  30  presents  open- loop  pole-zero  locations  for  fixed  Pq  = 3 deg 
and  trim  a<,  of  20,  and  21  deg*  Ihe  lateral  phugoid,  ojg|^,  steadily  decreases 
in  danqiing  as  a is  increased  and  moves  into  the  rij^t  half  plane  at  approximately 
21  deg  a.  The  dutch  roll  frequency  steadily  dv creases  and  damping  increases. 

The  short  period  is  unaffected  by  chsaige  in  a because  Ma  does  not  vary  in  this 
specific  a range  and  Zq,  was  not  varied  with  a in  this  simulation.  The  migra- 
tion of  poles  in  Fig.  30  is  generally  representative  of  all  nine  test  configurations. 

The  mid- frequency  zeros  of  9/8e  are  seen  to  consist  of  a complex  pair  which 
are  associated  with,  and  "track,"  the  dutch  roll  poles  as  a is  increased.  This 
again  is  typical  of  all  configurations  and  will  not  be  discussed  farther.  The 
other  zeros  of  the  pitch  numerator  include  a first  order  at  about  O.5  rad/sec 
which  is  relatively  unaffected  by  changes  in  a,  and  two  first- orders,  one  in 
the  right  half  plane  (1/T95),  which  drive  together,  become  complex,  and  move 
into  the  right  half  plane  at  approximately  20  deg  a. 

The  roll  numerator  initially  has  two  pairs  of  complex  zeros.  One  is 
associated  with  the  longitudinal  short-period  mode  and  essentially  ceincels 
that  pole.  The  other  is  the  conventional  cu(p.  As  angle  of  attack  is  Increased 
this  conqplex  zero  rapidly  moves  to  the  real  axis  where  it  splits  into  two  real 
zeros,  one  of  which  moves  into  the  ri^t  half  plane  at  approximately  20  deg  a. 

It  should  be  recalled  that  o\p  is  the  parameter  sometimes  referred  to  as  the 
aileron  control  divergence  parameter  and  is  considered  to  indicate  a tendency 
to  lateral  divergence  when  it  becomes  negative  (i.e.,  one  root  in  the  ri^t 
half  plane). 


*The  -rlean  aircraft  stalls  at  a = 21  deg  and  = 0 deg. 
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Figure  30,  Configuration  1 Open-Loop  Pole-Zero 
Migration  with  aj  Po  = 5 deg 


The  zeros  of  the  yaw  numerator  consist  of  one  real  and  two  complex  pairs. 
Again,  one  complex  pair  derives  from  the  coupled  longitudinal  short  period 
and  essentially  cancels  the  pole  of  this  mode.  The  other  complex  pair  is  the 
conventional  ojj.  which  lies  near  the  dutch  roll  pole  at  19  deg  angle  of  attack 
and  is  relatively  unaffected  by  the  2 deg  change  in  a shown  here.  The  real 
zero  is  the  conventional  l/Tj.  which  normally  is  of  approximately  the  same  vedue 
as  the  roll  subsidence  mode,  1/Tr. 

The  box  in  the  center  of  the  figure  presents  the  values  of  Npjjyjj  over  this 
a range  and  shows  that  this  parameter  remains  positive  until  even  higher  a is 
reached.  A three-dimensional  plot  of  versus  a and  p is  shown  in  Fig.  31. 

This  shows  that  if  p = 0 then  changes  sign  at  approximately  21  deg  a. 

However,  the  presence  of  a few  degrees  of  sideslip  will  keep  positive 

until  well  past  stall  a. 

At  this  point  the  piloting  technique  generally  employed  should  be  discussed. 
First,  in  attempting  to  track  and  stay  with  the  target  as  it  rami>ed  up,  pitch 
attitude  was  the  primary  loop  closure.  As  nose  slice  developed,  the  pilots 
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Figure  31.  Three-Dimensional  Plot  — - N^dyn  vs.  a,  p 
( Configuration  1 ) 

would  either  release  clLI  contiols  or  tend  to  hold  aileron,  rudder,  or  both 
for  several  seconds  at  a time  against  the  motion  of  the  Initial  slice.  Thus, 
the  lateral  controls  were  used  more  as  gross  moment  producers  to  oppose  air- 
craft motion  and  did  not  reflect  closed-loop  activity  In  the  usual  sense. 
When  aileron  or  rudder  was  used  in  eui  attempt  to  speed  reco\’ery,  pilot  com- 
mentary indicated  a general  problem  In  jdiaslng  either  control  to  tnprove 
recovery.  Generally,  the  use  of  these  controls  Increased  sideslip  and 
degraded  recovery-  For  exanqple.  Fig.  32  shows  a ”docile”  recovery  with  a 
complete  lack  of  aileron  or  rudder  ewtivity  in  both  departure  and  recovery. 
Figure  33&  indicates  how  use  of  rudder  and  aileron  can  aggravate  recovery 
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i and  also  the  extended  time  periods  full  contx^  deflection  is  used.  Three 

departure  and  two  recoveries  are  shown  in  these  traces.  Departure  was 
not  recoverable.  In  Departure  0 attempted  control  by  rudder  greatly 
excited  sideslip  (note  p and  r scale  changes  between  Figs.  ^2  and  3^)  but 
did  not  change  the  maximum  a that  could  be  achieved  before  recovery  bad  to 
be  initiated.  In  Departure  0 the  rudder  was  rapidly  bottomed  and  held. 

Recovery  was  achieved  with  stick  push  (to  a = 13  deg)  and  full  rudder.  In 
Departure  0 the  pilot  attempted  to  track  to  bi^^r  a and  possibly  was  a 
little  late  in  applying  corrective  rudder.  A small  amount  of  ailercn  was 
also  applied.  At  the  time  controls  were  released  the  sideslip  exceeded 
13  deg,  yaw  rate  was  greater  than  60  rad/sec,  and  a pitch  up  (inertia  coup- 
ling) was  underway.  Figure  33b  presents  pitch  attitude  and  heading  traces 
for  these  same  three  dei»rtures.  These  indicate  that  the  successful  recoveries 
had  peak  heading  excursions  of  approximately  3 unsuccessful  recovery 

was  initiated  when  heading  had  chan^d  about  20  deg.  One  second  later  the  yaw 

s was  30  deg. 

Recovery  with  all  controls  released  was  ;»ated  a 4 (Cooper-Harper)  by 
Pilot  JS  and  3 by  Pilot  ML.  Use  of  aileron  or  rudder  increased  the  average 
rating  to  7 for  JS  and  9 for  ML. 

f 

I Pilot  commentary  and  the  strip  charts  indicate  a "moderate"  nose  slice 

f 

i starts  at  a <.  21  deg  and  6 between  1 and  3 deg*  From  Fig.  50  there  are  two 

possible  parameters  associated  with  the  nose  slice.  One  is  the  longitudinal  4 

complex  zero,  tug.  The  other  is  the  roll  conplex  zero,  cn^j.  Both  result  in  , 

right  half  plane  zeros  at  approximately  20  deg  a.  On  the  basis  of  likely  loop 
closures  it  would  appear  that  the  pitch  zero,  cd0,  is  the  controlling  parameter. 

With  the  pitch-to-elevator  loop  closed  sufficiently  tight  to  control  the  short- 
period  attitude  tue  lateral-pliugoii  complex  pole  is  driven  into  tiie  complex 
zero  whicn  moves  into  the  right  half  plane  at  a 20  deg.  Thus,  the  departure 
s’  ouli  be  expectei  to  start  at  t!iat  angle  of  attack. 

As  a further  aid  in  understanding  the  mechanics  of  the  nose  slice,  closed- 
loop  modal  response  vectors  for  the  three  dynamic  modes  ( lateral-phugoid,  ugpj 
dutch  roll,  and  short  period,  ugp)  are  shown  in  Fig.  for  unity  gain 
pitch  attitude  clos'jre  and  21  deg  a.  This  shows  that  the  divergent  ugp  mode 
(the  closed-loop  pole  that  is  driven  into  the  right  half  plane  complex  zero,  ay.) 


Closed  Loop  Modal  Response  Vectors 


1 

* 


is  primarily  a roll-yaw  motion  (•?/’:  = I.U5)  in  whicl;  yaw  leads  roll  by 
approximately  30  deg  (positive  rotation  is  counterclockwise).  Althou^  the 
response  ratio  indicates  roll  is  greater  than  yaw,  in  the  simulation  heading 
displacement  should  be  more  apparent  to  the  pilot  than  bank  angle  because  a 
change  in  heading  of,  say,  3 deg  is  approximately  e quail  to  a semi- spam  of  the 
target  aircraft.  With  roll  lagging  yaw,  a bank  angle  deviation  of  5 deg 
would  not  be  as  apparent  because  of  the  lack  of  a direct  reference.  Thus, 
the  closed-loop  modal  response  coefficients  further  support  o)g  as  the  para- 
meter domlmuit  in  setting  nose  slice  chau'acteristics. 

Since  the  preferred  recovery  is  completely  open- loop,  modal  response 
vectors  for  the  19  deg  a trim  condition  au:e  presented  in  Pig.  35*  If  recovery 
is  initiated  (controls  released)  before  laurge  euaguleur  rates  develop,  inertia 
coupling  will  be  small  and  the  aircraft  will  return  to  this  trim  condition 
with  small  sideslip.  The  plots  show  that  the  laterad-phugoid,  is  again 

mainly  roll  and  yaw  (|cp/ilr|  = 5*  liU,  jcp/3l  = 22.3)  in  which  the  roll  Isigs  yaw 
by  approximately  lUO  deg.  Thus,  open  loop,  this  mode  has  considerably  greater 
rolling  motion  thau  closed  loopj  and  the  lag  between  and  cp  might  create  the 
control  phasing  problem  referred  to  1:^  the  pilots.  With  these  two  motions 
IUO  deg  out  of  phase  the  motion  would  be  characterized  by,  say,  left  yaw  and 
right  bank.  This  certainly  would  lead  to  difficulty  in  arresting  both  cp  and  \|r 
with  control  deflections  which  do  not  further  excite  sideslip  (i.e.,  without 
cross-control) . 

The  open- loop  dutch  roll  mode,  o)^j,  is  mainly  roll  and  sideslip  with 
|rp/p|  = 3.37  but,  interestingly,  Iq/p|  = 0.236  is  greater  than  [r/p|  = O.O83. 
Thus,  there  is  more  pitch  rate  than  yaw  rate  in  the  dutch  roll  mode.  The 
longitudinal  short  period,  cogp,  on  the  other  hand,  is  primarily  a pitch  euid 
roll  motion  with  cp/p  = 0.88  which  is  approximately  equal  to  the  0/a  ratio. 
Also,  p/a  = 1.71.  Thus,  the  greatest  longitudinal-lateral  coupling  is  in  the 
open- loop  short-period  mode. 

In  summary,  the  a at  which  Configuration  1 nose  slice  occurs  is  consistent 
with  movement  of  the  pitch  numerator  zero,  cog,  into  the  right  half  plane.  The 
closed-loop  modal  response  ratio  evaluated  at  this  zero  indicates  the  motion 
to  be  primarily  roll-yaw  in  which  the  yaw  would  be  more  apparent  in  a target 
tracking  task.  The  parameter  could  also  be  associated  with  nose  slice 
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initiated  at  or  slightly  greater  than  20  deg  a providing  the  sideslip  reaained 
very  near  zero.  However,  is  not  consistent  with  inst^ility  lAien  fi>  3 deg 

and  a < 22  deg  as  was  encountered  in  this  sianilation.  Ihe  preferred  recoveiy 
was  achieved  hy  release  of  all  controls  up  to  a = 21 .3  deg  and  p = 3 deg,  pro> 
viding  angular  rates  were  not  so  hl£^  at  release  as  to  induce  inertia  co\;9ling 
into  pitch.  T!ie  c^n>loop  modal  response  ratio  of  the  lateral-phugoid  is  con- 
sistent with  and  would  tend  to  predict  the  leasing  problems  encountered  in  the 
use  of  the  aileron  and  rudder  during  recovery.  If  the  pilot  attempted  to 
maintain  the  pitch  attitude  at  that  angle  at  which  nose  slice  started,  lateral- 
directional  control  could  not  be  maintained  and  severe  departure  resulted. 


2.  Confliuration  2 

Figure  36  presents  the  open-loop  pole>zero  survey  for  Configuration  2.  In 
this  configuration  I/T95  is,  initially,  deeper  into  the  riglit  half  plane  and 
is  unity.  same  range  of  values  as  in  Configuration  1 . 

The  values  of  remain  positive  and  in  the  left  half  plane  throughout  an 
angle  of  attack  range  exceeding  21  deg.  Movement  of  the  lateral-phugoid  is 
relatively  slight  between  19  and  20  deg  a but  rapidly  moves  into  the  ri^t 
half  plane  by  21  deg  a.  The  two  lower- frequency  zeros  of  the  e/6e  numerator 


Figure  36.  Configuration  2 Open-Loop  Pole-Zero 
Migration;  Pq  = 5 deg 
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come  together  in  the  rl^t  half  plane  and  the  resulting  cooplex  zeros 
remain  in  the  ri^t  half  plane  at  all  times. 

The  pilot  cooraented  that  this  configuration  dei>arts  quickly  starting 
at  a < 21  deg  and  the  slice  speeds  as  builds.  This  configuration  also 
would  recover  if  the  controls  were  released  prior  to  a = 21 .5  deg  and  p 
2 to  5 deg.  It  did  not  recover  if  the  sideslip  was  allowed  to  build  to 
the  order  of  5 deg.  Problems  were  again  encountered  in  coordinating  use 
of  aileron  or  rudder  and,  in  general,  they  aggravated  recovery.  Recovery 
with  release  of  all  controls  yielded  a rating  of  by  JS  and  7 by  ML. 

It  may  be  noted  from  Fig.  36  that  both  and  Np^yn  remain  positive  at 
the  a and  P values  where  departure  was  obtained.  Again,  only  the  pitch 
parameter,  cog,  correlates  with  the  departure  conlitions. 

3 . Configuration  3 

For  this  configuration  the  lateral-phugoid  was  changed  in.o  the  more 
conventional  spiral  and  roll  subsidence  modes  with  1/Tr  = 0.46j  Np^jyjj  was 
increased  to  1 .58;  1 /Te^  remained  at  approximately  the  same  value  as  in 
Configuration  2j  and  av,'Au(i  was  returned  to  the  base  (Configuration  1 ) A-7 
value. 

The  open-loop  pole-zero  migrations  with  a are  shown  in  Fig.  37.  The 
major  difference  between  this  and  the  previous  case  is  that  the  lateral- 
phugoid,  which  results  from  a rejoining  of  the  spiral  and  roll  subsidence 
at  about  20  deg  a,  remains  well  in  the  left  half  plane  at  21  deg  a. 

Pilot  commentary  indicated  this  configuration  exhibited  a "mild"  slice 
characteristic  which  could  be  controlled  to  approximately  20  deg  a but 
departed  for  a > 20  deg.  Recovery  could  be  achieved  by  releasing  all  con- 
trols iip  to  22  deg  a and  b deg  p and,  with  stick  push  forward  of  neutral,  to 
20  deg  a and  13  deg  p.  Use  of  lateral  control  would  rapidly  build  sideslipj 
and  this  configuration  generally  exhibited  larger  sideslips  than  the  previous 
case.  Both  pilots  rated  the  departure /recovery  a 3 for  release  of  all  con- 
trols. Stick  push  recovery  after  development  of  large  sideslip  was  rated  a 7. 

Again  for  Configuration  T both  lateral  parameters  nir.  and  remained 

positive  (with  comfortable  margin)  and  the  open-loop  lateral-phugoid  remained 
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Figure  37.  Configuration  5 Open-Loop  Pole-Zero 
Migration;  Pq  = 5 deg 


stable  at  the  angle  of  atx;ack  at  which  nose  slice  was  exhibited.  Only  cug 
lies  in  the  right  hsdf  plane  and  provides  a good  indication  of  closed-loop 
divergence  tendencies. 

4.  Configuration  k 

This  configuration  was  intended  to  provide  docile  handling  at  high  angle 
of  attack.  To  accomplish  this  ’’le  open-loop  parameter  l/Tg^  was  moved  to  the 
left  half  plane  and  uip/o^  was  set  to  unity.  The  survey  plots  are  presented 
in  Fig,  38.  As  with  the  basic  (Configuration  1 ) A-7,  the  lateral-phugoid 
steadily  moves  toward  and  into  the  right  half  plane  at  approximately  21  deg  a. 
The  two  low-frequency,  real  zeros  of  the  pitch  numerator  couple  at  between 
19-20  deg  a Md  again  move  into  the  ri§^t  heuLf  plane  at  between  20-21  deg  a. 
The  aileron  zero,  ouip,  remains  in  the  left  half  plane  said  Np,j|yn  remains  posi- 
tive throughout  this  a range. 

Pilot  cot,imentary  indicated  a "slow”  rate  of  nose  slice  was  obtained  at 
about  21  deg  a but  in  this  instance  the  pilot  could  hold  up  to  23  deg  a for 
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Figure  58.  Configuration  k Open- Loop  Pole-Zero 
Migration;  Po  = 5 deg 


several  seconds  while  using  lateral  control  to  contain  the  divergent  nose 
wander  within  acceptable  bounds  (see  Fig.  59) • This  often  could  be  accom- 
plished without  exciting  "large"  sideslip  and  recovery  could  be  achieved  by 
merely  releasing  all  control.  This  configuration  was  considered  very  docile 
(as  planned)  unless  the  pilot  induced  large  sideslip  while  attempting  to 
maintain  control  at  the  high  angle  of  attack.  Again,  both  pilots  rated  this 
a Cooper-Harper  of  3 for  recovery  involving  only  release  of  controls.  If 
lateral  control  were  en^jloyed  following  longitudinal  stick  release,  problems 
were  again  encountere  i in  phasing  the  lateral  control  to  iirprove  recovery  and 
ratings  suffered  (average  of  7 to  8). 

The  nose  slice  again  correlates  with  u)g  moving  into  the  right  half  plane. 

It  is  particulsurly  interesting  though  that  this  configuration  could  be  flown 
for  several  seconds  at  angles  of  attack  at  or  above  22  deg  a since  this  should 
result  in  a closed-loop  instability  as  in  the  previous  cases,  especially  since 
%dyn  approach  negative  values.  An  explanation  for  this  relatively 
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docile,  Moi-controllable  brtiavlor  at  hl^er  angles  of  attack  is  that  the 
initial  magnitude  and  slope  (with  p)  of  Nc  are  reduced  hy  s factor  of  ^ 

(from  -1 .6  to  -0.5  at  p * 5 deg)  in  order  to  move  l/Tej  into  the  left  half 
plane.  Thus,  there  is  a significant  i*eduction  in  longj.tudinal  to  direc- 
tional coupling.  This  then  requires  less  lateral  control  activity  and  hence 
produces  less  sideslip  excitation.  Departure  is  delayed  hut  not  prevented. 

; . Cooflfuratlon  ^ 

The  survey  plot  for  Configuration  5 is  presented  in  Fig.  UO.  Again,  l/Tej 
is  in  the  rig^t  half  plane.  The  principal  change  from  Configuration  4 is  an 
increase  in  which,  among  other  effects,  increases  Aileron  yaw, 

was  also  returned  to  the  nominal  (Configuration  1 ) A-7  value}  and  these  two 
chiuiges  resulted  in  ui(p/a^  = 0.534,  which  shouid  prc/ide  roll  reversal  tenden- 
cies. The  increase  in/p  results  in  considerably  less  migration  of  all  poles 
and  zeros  (except  ui^p)  for  the  2 deg  a change  shown  here.  The  combination  of 
positive  1 /T05  and  increased /p  resulted  in  the  most  stable  configuration  of 
those  investigated  (no  poles  or  zeros  in  the  right  half  plane  at  stall)  and 


Figure  40.  Configuration  5 Open-Loop  Pole-Zero  Migration 


also  provided  the  least  change  In  vehicle  dynamics  as  stall  Is  approached. 

Both  (ogjt  and  <hq  reflect  a stable  lateral-phugold  at  stall  (a  = 21  deg). 

Althou^  transfer  function  roots  were  not  calculated  for  a's  greater  than 
stall,  the  a at  which  oig  reaches  the  rigiit  half  plane  may  be  estimated  by 
comparing  the  root  movements  for  Configurations  U,  and  6.  The  only 
significant  difference  in  these  configurations  Is  a change  InjC^. 

Configuration  6:  jCp  = —1.75  , Oojg  = 19.5-2O  deg  (Fig.  U2) 

Configuration  U;  jCp  = -^  , = 20.5-21  deg  (Fig.  38) 

Configuration  5:  = -5.55 

By  compturison,  it  would  be  expected  Ocug  (the  a at  which  odq  crosses  into  ri^t 
half  plane)  for  Configuration  5 to  be  about  21 .5  to  22  deg. 

Pilot  coomentary  indicates  this  configuration  is  still  under  control  at 
21  deg  a and  slow  nose  slice  is  not  encountered  \intil  approximately  22.5  to 
23  deg  a.  It  could  be  taken  to  23  deg  a and  6 deg  p or  more  and  still 
recovered  with  only  release  of  controls.  One  pilot  rated  this  a the 
other  pilot  rated  it  a 3.  There  were  completely  different  opinions  regard- 
ing the  effectiveness  and  usefulness  of  aileron  and  rudder.  One  pilot  thou^t 
the  rudder  effective  while  aileron  was  not;  the  other  had  the  reverse  opinion. 
However,  both  could  maintain  control  quite  well  while  at  the  elevated  angles 
of  attack.  One  pilot  commented  that  he  had  so  much  confidence  at  the  hi^er 
a's  that  he  would  get  into  trouble  by  attempting  to  track  too  far  into  the 
stall  region  (see  Fig.  4l).  He  also  considered  the  departure  warning  to  be 
inadequate.  It  was  difficult  for  either  pilot  to  decide  when  to  "give  up" 
and  attempt  recovery  because  lateral-directional  control  remained  quite 
effective.  However,  sideslip  would  often  continue  to  build  up  until  finally 
a severe  departure  would  result  (e.g.,  the  second  stall  in  Fig.  4l).  Because 
there  was  a tendency  to  go  too  far,  both  pilots  downgraded  the  departure/ 
recovery  rating. 

In  summary,  this  configuration  proved  to  be  the  most  stable  and  difficult 
to  d^>art,  as  planned,  but  proved  surprising  in  that  the  warning  was  considered 
inadequate.  Ibe  Importamt  aspect  remains,  however,  that  the  presence  (or  lack 
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of)  nose  slice  at  a given  a Is  consistent  with  the  location  of  the  ccnqplex 
zero  (U0  cukL  that  this  configuration  Which  vas  selected  to  provide  the  best 
hl£^  a characteristics  did  Indeed  prove  the  most  difficult  to  obtain 
departure. 

6-  Conflguratloa  6 and  7 

The  survey  plots  representing  both  Configurations  6 and  7 are  presented 
in  Fig.  42.  The  change  from  Configuration  5 to  Configuration  6 is  accom- 
plished by  changing  one  derivative,  from  -5.55  to  —1 .75.  Configuration  7 
differs  from  Configuration  6 in  that  the  coupling  derivative, is  reduced 
from  1.55  to  0.75  for  Configuration  7 to  determine  whether  this  influences 
the  roll  characteristics  of  the  departure.  It  did  not  produce  a significant 
difference  in  departure/recovery  characteristics  or  the  pole-zero  locations 
shown  in  Fig.  42.  These  initially  correspond  to  a stable  l/Tej,  o^/a^  = 0.63, 
low-frequency  lateral  characteristics  separated  into  a spiral  and  roll 
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Figure  42,  Configurations  6 and  7 Open-Loop  Pole-Zero  Migration 
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subsidence,  and  an  initial  Np^jyn  = O.56  iidiich  becomes  negative  at  approxi- 
mately 20.5  deg  a.  Increasing  «ingle  of  attack  results  in  relatively  large 
pole-zero  migrations  from  these  initial  values.  The  roll  subsidence/spiral 
rapidly  couple,  pass  into  the  ri^t  half  plane  at  20  deg  a,  and  produce  two 
first-order  divergence  modes  (RHP)  by  21  deg  a.  The  dutch  roll  mode  damping 
increases,  becomes  critically  dattqped  at  approximately  21  deg  a,  and  separates 
into  two  first-order  stable  poJes.  The  low-frequency  zeros  of  track 
almost  precisely  with  the  roll  subsidence/spiral  poles  and  thus  produce  a 
complex  ri^t  half  plane  pair  by  20  deg  a.  In  the  roll  numerator,  moves 
to  the  real  axis  and  produces  a first-order  ri^t  half  plane  zero  at  approxi- 
mately 19.5  deg  a and  this  zero  moves  further  into  the  ri^t  half  plane  as  a 
increases.  Only  the  yaw  zeros  remain  in  the  left  half  plane. 

This  configuration  produced  a very  flat-appearing  yaw  departure  at  a 
greater  than  20  deg,  accompanied  by  a very  large  and  rapid  buildup  of  sideslip. 
One  recovery  was  obtained  from  21.5  deg  a and  h.'J  deg  pj  however,  recovery  was, 
in  general,  not  achieved  unless  back  stick  was  released  upon  first  detection 
of  nose  slice.  Following  release  of  back  stick,  one  pilot  thought  aileron 
improved  recovery.  Use  of  rudder  in  recovery  was  considered  difficult.  One 
pilot  rated  control-free  recovery  a 4 while  the  other  rated  it  between  a 7 
6«id  9-  The  major  difference  appeared  to  lie  in  how  rapidly  the  pilot  initiated 
recovery  once  the  nose  slice  started. 

The  more  positive  and  rapid  departizre  obtained  with  this  configuration  was 
apparently  related  to  the  almost  simultaneous  achievement  of  negative  values 
with  the  three  parameters  co0,  cu^,  and 

7.  Configuration  6 

The  survey  plots  for  Configuration  8 are  shown  in  Fig.  1+3 . In  this  case 
an  initial  condition  of  Np^jy^  = 0 was  achieved  by  retaining  the  low  Xp  (Con- 
figurations 6 and  7)  and  setting  Np  = -0.5.  This  results  in  an  initial 
trinmed  condition  where  the  lateral -phugoid  mode  is  unstable.  Further 
increase  in  a results  in  two  first-order  divergences,  one  of  which  is  very 
rapid.  Both  the  pitch  and  roll  numerators  also  start  with  zeros  in  the  right 
half  plane)  thus,  neither  loop  closure  will  stabilize  the  divergence.  As 
would  be  expected,  this  configuration  produced  immediate  divergence  which 
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Figure  43.  Configuration  8 Open-Loop  Pole-Zero  Migration 


could  not  be  arrested  by  use  of  control  unless  longitudinal  stick  was  pushed 
Md  held  forward  of  neutral  until  recovery.  This  reduced  the  angle  of  attack 
below  19  deg  *nd  returned  the  lateral-phugoid  mode  to  left-half -plane  sta- 
bility. Very  little  could  be  done  with  this  configuration,  and  recovery  had 
to  be  initiated  upon  the  first  indication  of  the  nose  slice.  Ihe  use  of 
lateral-directional  control  made  this  configuration  essentially  unrecoverable. 
While  it  could  be  recovered  maintedning  forward  longitudinal  stick,  this  did 
require  more  effort  and  concentration  on  the  part  of  the  pilot  and  therefore 
degraded  the  recovery  rating  to  a 6 or  7 for  both  pilots. 


8.  Configuration  9 


The  final  configuration,  shown  in  Fig.  44,  offered  improved  initial  low- 
frequency  lateral  dynamics  (l/5^  = 0.53j  1/Tbq  = 0.03)  but  decreased  dutch 
roll  (5^  = 0.034)  over  that  of  the  basic  A-7.  IMs  was  accooqjlished 

by  increasing  Np  and  making  the  term  (Kp  — g/Uo)  negative  (i.e.,  decreasing 
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Figure  44.  Coofiguratir-n  ) Open-Loop  Pole-Zero  Migration 

the  total  yaw  due  to  roll  rate).  is  returned  to  its  initial  (Configura- 
tion l)  value j therefore,  1 /T0j  is  negative. 

The  princii>al  effect  of  an  open- loop  increase  in  a is  to  cause  the  dutch 
roil  to  diverge.  Movenent  of*  all  other  poles  is  relatively  small.  !Hie  low- 
frequency  zeros  of  N|g  couple  into  the  complex  009  which  remains  in  the  left 
half  plane  at  21  deg  a but  would  appear  likely  to  move  to  the  ri^t  half  plane 
shortly  thereafter.  Ihe  complex  zeros,  0)^^,  remain  in  the  left  half  plane  and 
move  toward  the  real  axis.  In  pitch  und  yaw,  the  dutch  roll  ccMnplex  pole-zero 
pairs  essentially  cancel  so  there  is  little  evidence  of  the  dutch  roll  in 
these  motions.  The  dutch  roll  is  therefore  primarily  a roll  motion.  If  the 
pilot  attempts  to  control  or  dan^)  dutch  roll  wxth  ailerons  at  the  hi^er  a's 
we  would  expect  the  closed-loop  root  migration  to  be  as  shown  in  the  following 
sketch.  Here  the  roll  subsidence  and  spiral  drive  into  u)(p  while  u:^  moves  to 
higher  frequency,  ihe  closed-loop  dutch  roll  also  is  further  destabilized  by 
the  pilot's  time  delay  and  other  hi^-frequency  lags.  !Hius,  this  configuration 
mi£^t  be  expected  to  exhibit  a wing  rock  tendency  prior  to  nose  slice. 
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Pilot  connentaxy  and  strip  chart  records  Indicate  nose  slice  did  steort 
at  about  22  deg  a.  One  pilot  described  it  as  having  a roll-off  tendency 
before  the  departure  angle  of  attack  is  reached.  The  other  described  it  as 
a tendency  to  wing  rock.  Recovery  was  achieved  by  release  of  all  controls 
at  p tq)  to  ±7  deg  and  a = 22  deg.  Use  of  siny  lateral  control  only  aggravated 
the  recovery.  This  configuration  was  rated  a 3 for  inmediate  release  of  all 
controls.  If  an  attaiqpt  was  made  to  damp  the  dutch  roll  motion  during 
departure  and  recovery  it  was  rated  between  6 and  10,  depending  upon  how 
aggressively  the  pilot  opposed  the  roll  motion. 

9.  Bmmry  Md  Cenoluiloni 

It  is  concluded  that  the  initial  and  final  nose  slice  characteristics 
and  pilot  comments  best  relate  to  the  i«rameters  1/Tqj  or  o)©.  The  angle 
of  attack  at  which  departure  occurs  can  be  increased  significantly  by  elimina- 
tion of  N^,.  The  conventional  latereLL  parameters,  Np,jyjj  and  oj^),  certainly  con- 
tribute to,  but  are  not  the  key  factors  in,  the  A-7  nose  slice.  It  appears 
that  when  these  x>arameters  become  negative  at  the  same  time  the  longitudinal 
parameters  are  in  the  ri^t  half  plane  the  rapidity  and/or  severity  of  nose 
slice  dep€u*ture  is  increased.  It  is  also  apparent  that  lateral-directional 
control  is  difficult  to  coordinate  during  the  departure  and/or  recovery  phases 
and  generally  aggravates  rather  than  enhances  recovery.  If  system  dynamics 
and/or  control  characteristics  do  not  change  appreciably  as  a increases  the 
masked  departure  onset  can  result  in  attempts  to  maintain  control  to  higher 
a and,  eventually,  a more  severe  departure.  One  key  to  successful  recovery 
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with  minijnal  pilot  workload  or  compensation  is  that  the  aircraft  be  trimmed, 
stick  force  zero,  at  an  angle-of«attack  for  >ihich  the  lateral-phugoid  and 
dutch  roll  are  stable.  1iierefo3*e,  the  open»loop  parameter  is  a key 

factor  in  recovery.  Finally,  the  results  siqpport  the  concept  that  wing  rock 
is  related  to  an  unstable  dutch  roll  which  is  farther  aggravated  by  pilot 
attenpts  at  roll  control. 

C.  EFFECT  OF  PILOT  TECHHIQUE 

The  pilot  ratings  assigned  to  the  different  configurations,  according  to 
the  Cooper-Harper  handling  quality  rating  scale  of  Pig.  45,  varied  widely 
with  pilot  technique  (aigressiveness,  anti-  versus  pro-departure  recovery 
control,  stick  release  versus  stick  push  recovery,  etc.).  In  general  any 
attempt  at  use  of  lateral-directional  control  to  reduce  departure  excursions 
resulted  in  a significant  increase  in  numerical  rating.  This  was  traced  to 
several  factors.  First,  there  was  a tendency  to  use  less  lateral-airectional 
control  in  opposition  to  the  initial  departure  motion  if  recovery  was  to  be 
achieved  via  release  of  all  controls.  Conversely,  if  the  task  called  for 
use  of  aileron  or  rudder  during  recovery,  there  was  a strong  tendency  to 
significant  use  of  this  control  in  the  pre-departure  tracking  task.  Effec- 
tive application  of  such  control  required  significant  pilot  concentration 
(compensation)  and  hence  degraded  the  rating.  The  application  of  latered 
control  to  oppose  the  initial  departure  motion  either  allowed  deeper  stall 
penetration  or  increased  sideslip  — or  both.  This  then  resulted  in  a more 
violent  initial  nose  slice  and,  quite  often,  a "secondary"  departure  due  to 
roll- yaw  inertia  coupling  into  pitch. 

As  indicated  previously,  the  pilots  were  requested  to  assess  or  evaluate 
two  separate  types  of  recovery.  One  might  be  called  aggravated  control 
recovery  in  which  an  attempt  was  made  to  hold  the  departure  pitch  attitude 
and  maintain  lateral-directional  control  of  the  aircraft.  This  was  not 
sucessful  in  any  of  the  cases  and  always  resulted  in  a rating  of  10,  The 
second  method  was  to  release  and/or  push  forward  on  the  stick  to  reduce 
angle- of- attack.  The  latter  technique  was  further  categorized  by  specific 
use  of  lateral  control  following  release  of  back  stick.  The  pilot  ratings 
obtained  for  the  matrix  of  dynamic  configurations  and  recovery  techniques 
are  given  in  Table  B:  recovery  A identifies  release  of  longitudinal  stick 
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TABLE  8.  SUMMARY  OF  PILOT  RATINGS  FOR  RECOVERY 
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only  with  no  lateral-directional  inputs  following  stick  release,  B is  use 
of  aileron  only,  C is  rudder  only,  D is  aileron  and  rudder,  and  E is  rudder 
only  but  withjCg^  reduced  by  a factor  of  four.  The  latter  resulted  in  very 
nearly  pure  yaw  rudder  for  all  configurations. 

In  addition  to  these  scheduled  variations  in  control  technique  there 
were  several  unplanned.  One  pilot  was  more  aggressive  in  attoopting  to 
track  the  target  to  high  a's  before  giving  up  and  Initiating  recovery. 

(The  pilot  measurements  reported  in  ihe  following  sections  showed  this 
pilot  to  be  operating  at  the  higher  gain  of  the  two.)  This  pilot  also 
routinely  used  forward  stick  in  evaluating  various  recovery  techniques. 

The  other  pilot  did  not  like  to  use  forward  stick  and  appeared  to  use  this 
recovery  technique  only  as  a last  resort.  However,  he  routinely  investigated 
application  of  c<xrtr61s  with  the  direction  of  nose  slice  as  well  as  control 
against.  Control  "with"  sometimes  proved  quite  successful  in  converting  a 
nose  slice  departure  into  a controlled  36O  deg  roll  maneuver. 

Table  8 indicates  the  wide  range  in  ratings  achieved  for  the  configurations 
and  recovery  techniques.  The  search  for  the  cause  of  the  wide  rating  dis- 
crepancies between  the  two  pilots  for  any  given  configuration  and  recovery 
technique  led  to  identification  of  the  individual  pilot  techniques  noted 
above.  For  exanq>le,  the  first  column  indicates  by  a Y the  instances  in 
which  it  could  be  identified  from  the  time  traces  that  Pilot  JS  successfully 
converted  a nose  slice  into  a 56O  deg  roll.  These  almost  always  resulted  in 
a 3 to  5 point  inprovement  in  rating.  There  were  other  attenpts  at  other 
flight  conditions  which  were  unsuccessful  and  are  not  identified  in  this 
column.  Those  instances  when  this  pilot  used  forward  stick  to  effect  recov- 
ery are  indicated  in  the  stick  push  column.  Those  cases  left  blank  correspond 
to  a "normal"  stick  release;  there  was  no  attempt  to  use  forward  stick  except 
as  a last  resort.  Pilot  ML  routinely  evaluated  both  a stick  release  recovery 
and  a stick  push  recovery,  therefore,  pilot  ratings  for  both  techniques  are 
given. 

In  essentially  every  instance  both  pilots  rated  recovery  technique  A 
superior  to  any  other.  There  appesured  to  be  several  factors  involved  here. 
First,  it  is  the  simplest  technique  and  requires  no  compensation  or  workload 
on  the  part  of  the  pilot.  Second,  as  indicated  previously,  there  was  a 
tendency  for  less  use  of  aileron  or  rudder  to  oppose  the  initial  departure 
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motion  (or  Just  give  up  sooner)  and  hence  the  post- release  aircraft  motlcm 
was  less  violent,  there  was  less  chance  for  significant  inertia  coupling 
gyrations,  and  positive  dall^>ing  of  aircraft  pitch  wd  roll  motion  with 
respect  to  the  horizon  was  more  readily  aiq^ent.  This  approach  to  dq>arture 
thus  takes  the  appearance  of  a deliberate,  straight-ahead,  1 g stall  and  the 
ratings  obtained  may  be  more  applicable  to  such  a maneuver  rather  than  ah 
inadvertent  departure  from  tracking. 

For  other  departure  characteristics  and  recovery  techniques  (i.e,  B, 

C,  D,  E, ) the  ratings  are  highly  variable  with  so  many  interacting  factors 
that  few  trends  can  be  determined.  However,  the  generally  more  aggressive 
application  of  control  in  the  early  departure  sequence  may  make  these  ratings 
more  appropriate  to  actual  inadvertent  departure  and  recovery.  Attenpts  to 
correlate  ratings  with  individual  parameter  variations  proved  unsuccessful. 

It  was  possible  to  resolve  most  of  the  large  (4  to  5 point)  discrepancies 
between  the  two  pilots  on  the  basis  of  the  previously  mentioned  360  deg  roll 
and  stick  push  recovery  techniques.  The  Configuration  8 ratings  show  that 
an  aircraft  trimmed  hands-off  to  a dynamically  unstable  attitude  is  non- 
recoverable  unless  the  pilot  is  willing  to  use  full  forward  stick.  However, 
Configuration  5 which  was  basically  the  most  stable  and  docile  configuration 
also  proved  nearly  as  non- recoverable  for  Pilot  ML,  because  there  was  no  clear 
indication  of  control  loss  and  he  would  continue  tracking  the  target  until 
one  or  more  controls  bottomed.  Relaxing  controls  to  neutral  at  that  point 
gave  a violent  departure. 

The  pilot  ratings  obtained  were  inherently  influenced  by  the  combined 
departure  warning,  departure  (closed-loop)  characteristics,  and  recovery 
(open-loop)  properties.  These  could  not  be  adequately  separated,  and  in 
hindsight,  there  is  question  whether  the  Cooper-Harper  handling  quality  rating 
scale  is  really  applicable.  For  example,  the  Cocker- Harper  rating  requires  a 
first  assessment  as  to  whether  the  aircraft  is  controllable  or  not  controllable. 
If  controllable,  the  decision  is  whether  performance  is  adequate  or  not.  Depar- 
ture involves  an  inherent,  short  term  loss  of  control.  It  is  first  a matter 
of  recognition  and  then  a matter  of  aircraft  response  — Did  the  pilot  think 
he  still  had  control?  Did  he  indeed  have  seme  margin  of  control?  How  much 
aircraft  motion  was  involved  before  the  pilot  regained  control?  Did  the  pilot 
regain  control?  How  much  concentration  was  required  to  detect  and  take  action? 
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All  of  these  factors  or  conslderatioos  Influenced  each  of  the  single  pilot 
ratings  obtained  here.  However,  the  overall  gyration  Involves  a succession 
of  Interrelated  events  too  broad  In  scope  to  relate  to  a sln^e  set  of 
descriptors.  Thus  separate  ratings  and  descriptors  for  warning,  d^arture, 
and  recovery  appear  necessary. 

It  should  also  be  noted  that  all  configurations  enplqyed  here  exhibited 
good  pitch  control  and  response  characteristics.  The  aircraft  did  not 
exhibit  longitudinal  Instabilities  except  as  a result  of  Inertial  coupling. 

The  Initial  loss  of  control  always  was  In  the  lateral-directional  modes. 

Recovery  always  depended  upon  reducing  the  angle  of  attack.  A rapid  and 
positive  pitch  response  to  control  release  or  stick  push  was  central  to 
the  recoveries  obtained  and  the  pilot  ratings  obtained.  When  any  attenpt 
was  made  to  hold  pitch  attitude  or  angle  of  attack  constant  after  the 
Initial  nose  slice  was  detected,  control  of  lateral-directional  motion 
could  not  be  maintained  and  a violent  departure  resulted  from  which  recovery 
was  inpossible.  If  a configuration  were  enployed  which  exhibited  pitch  up 
or  loss  of  longitudinal  control  effectiveness  as  stall  was  approached,  the 
results  might  be  considerably  different  from  those  obtained  here. 

D.  SUMABX  AND  CONaUSICNS 

Nine  airframe  dynamic  configurations  were  investigated  to  assess  various 
high  angle  of  attack  departure  and  recovery  characteristics.  All  exhibited 
nose  slice  type  departure  and  Configuration  9 also  had  a wing  rock  tendency. 

In  all  cases  the  a at  which  nose  slice  initiated  or  became  apparent  to  the  pilot 
correlated  most  closely  with  the  pitch  attitude  numerator  conplex  zero,  cug, 
moving  into  the  ri^t  half  plane.  Once  the  departure  was  triggered,  safe 
recovery  was  achieved  if  the  aircraft  was  trimmed  hands  off  for  an  a at 
which  was  greater  than  zero,  and  all  controls  were  released  before; 

a.  Large  anti- departure  aileron  and/or  rudder  deflection 
had  been  applied 

b.  A’chieving  approximately  ±5  deg  P or  i|f 

c.  Achieving  approximately  ±10  deg/sec  r 

Some  configurations  were  recoverable  from  larger  excursions  of  P or  \|f,  however, 
recovery  could  not  always  be  obtained.  The  use  of  aileron  or  rudder  to  eppose 
nose  slice  almost  always  aggravated  the  departure,  and  the  use  of  these  controls 
during  recovery  often  decreased  the  chance  of  recovery. 
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Departure  onset  did  not  appear  correlated  with  (d,,  or  No  * but  these 

““  ^ ^dyn 

parameters  did  influence  departure  severity  and  recovery  characteristics. 

If  all  parameters  (1/T0  , cuq,  (%,  o)gn,  Ng  ) were  in  the  right  half  plane 

or  negative,  rapid  departure  was  obtained  at  low  a and  with  little  warning 

and  hence  little  time  for  the  pilot  to  react.  If  the  lateral-directional 

parameters  (o^,  oigp,  remained  stable  (left  half  plane  or  positive), 

recovery  was  more  positive  and  the  aircraft  could  penetrate  further  into 

stall/departure  and  be  recovered,  however,  departure  warning  was  less 

pronounced. 

The  pilots  apparently  used  the  change  in  aircraft  dynamics  with  increasing 
a and  P as  departure  warning.  The  configuration  having  the  best  lateral- 
directional  characteristics  but  also  the  least  changes  in  dynamics  with  a 
was  downgraded  severely  because  of  inadequate  warning.  In  this  respect 
a steady  degradation  in  handling  and  a slow  but  firm  nose  slice  provided 
the  best  warning  and  time  to  react. 

Data  and  pilot  ratings  obtained  for  the  stick  release  only  recoveries 
probably  relate  mostly  to  1 g,  straight-ahead  stall  maneuvers.  With  few 
exceptions  these  involved  minimum  penetration  into  stall/departure,  were 
reswiily  recoverable  control  free,  and  received  Cooper-Harper  ratings  of 
3 to  6 from  one  pilot  and  3 to  7 for  the  other  (using  stick  push).  Data 
and  ratings  obtained  for  all  other  types  of  recoveries  (e.g.,  B,  C,  D) 
probably  relate  better  to  actual  tracking  and  maneuvering  situations  in 
which  inadvertent  departure  — or  delay  of  departure  — is  encountered. 

These  cases  consistently  resulted  in  Cooper-Harper  ratings  of  6 to  10.  The 
very  significant  difference  between  these  two  sets  of  ratings  may  reflect 
actual  fli^t  situations  in  which  stall  and/or  spin  characteristics  obtained 
from  intentional  entries  (including  application  of  pro- spin  controls)  tend 
to  be  considerably  milder  than  those  obtained  from  inadvertent  entries. 

The  resulting  pilot  ratings  could  not  be  correlated  with  any  single  flying 
quality  parameter.  The  departure  chatracteri sties  were  most  dq>endent  upon 
closed  loop  (i.e.,  pitch  numerator)  parameters  with  some  influence  from  open 

*Ref.  18  also  cites  lack  of  correlation  of  nose  slice  and  Na , or  cu^ 

. , ^dyn  <P 

parameter  values. 
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loop  parameters.  Where  recovery  was  possible  it  was  most  dq>endent  upon 
airframe  open-loop  parameters  (including  trim)  but  the  possibility  of 
successful  recovery  was  highly  dependent  uixm  pilot  technique^  aggressiveness^ 
sensitivity  to  cues,  etc. 

All  configurations  investigated  here  exhibited  good  pitch  ccmtrol  and 
response  characteristics.  If  a configuratiai  was  enployed  which  ejdiibited 
pitch-up,  loss  of  longitudinal  control  effectiveness,  or  any  other  in^ir- 
ment  of  l<xigltudinal  control,  the  results  ml^t  be  cmsiderably  different 
from  those  obtained  here. 

Finally,  it  should  be  noted  that  the  conventional  closed-loop  pilot  model 
does  not  appear  appropriate  for  analysis  of  departure  and  recovery.  Immedi- 
ately prior  to  departure  the  lateral- directional  control  (if  employed  at  all) 
tended  toward  steady  or  trim  type  applicatim  against  the  motioi  of  the  nose 
slice.  Following  departure  the  preferred  method  of  recovery  was  to  release  all 
controls.  Thus  the  pilot  provided  no  inputs  and  the  recovery  was  achieved 
through  qpen-loop,  airframe- alone  stability  and  danping.  Again,  if  control 
was  applied  during  the  recovery,  the  preferred  technique  was  to  hold  constant 
foivard  stick  to  assure  decrease  in  a and  to  prevent  pitch -up  due  to  inertia 
coupling.  As  the  aircraft  motion  subsided  the  stick  was  returned  to  neutral. 
Thus  a closed  loop  pilot  model,  per  se,  does  not  appear  appropriate  for 
departure/recovery  application. 
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Multiloop,  multi-irq)ut,  multi-output,  pilot-vehicle  describing  function 
measurements  were  made  using  three  different  airframe  dynamic  configurations 
representative  of  high-angle-of-attack,  pre-stall',  flight  and  a control  task 
closely  approximating  air-to-air  tracking.  One  configuration  assumed  sym- 
metric flight  (3=0)  and  hence  had  uncoupled  lateral-longitudinal  dynamics. 

Two  reflected  coupled  lateral-longitudinal  dynamics  discussed  in  the  previous 
section  (Configurations  1 and  3).  ihree  control  loop  structures  were  en^loyed 
as  previously  shown  in  Fig.  27.  Two  were  two  loop  pitch-roll  tracking  tasks 
with  the  heading  degree  of  freedom  eliminated  to  prevent  inadvertent  yaw  error 
from  biasing  the  pilot's  control  of  roll  attitude.  As  indicated  in  Section  III, 
the  command  forcing  function  was  introduced  one  axis  at  a time  as  the  primary 
tracking  task.  The  secondary  task  was  then  related  to  the  inter-axis  coupling 
of  the  airframe.  Since  the  principal  coupling  was  from  the  longitudinal  into 
lateral  motion,  the  pitch  tracking  task  provided  the  highest  workload  for 
Configurations  1 and  5.  The  third  loop  structure  was  the  most  realistic  of 
the  three.  In  this  instance  the  forcing  function  commanded  b6Uik  of  the  target 
aircraft  which  was  also  free  to  change  heading  (and  pitch).  The  pilot  was 
instructed  to  track  heading  and  pitch  (i.e.,  to  fly  in- trail).  In  order  to 
test  the  concept  that  pilots  use  target  bank  angle  as  a primary  cue,  several 
heading  tracking  runs  were  made  with  target  wings  and  tail  removed  from  the 
display.  It  was  deteruiined  that  this  path  was  of  little  consequence  on  the 
describing  function  measures;  therefore,  the  pilot  was  apparently  making  little 
use  of  target  bank  angle  cue  in  this  particular  set  of  tracking  tasks. 

Describing  function  measurements  were  made  on  both  pilots.  Each  pilot 
was  edlowed  two  training  sessions  of  roughly  one-half  day  each  to  learn  to 
track  the  target  for  100  sec  without  inducing  a depeurture.  During  all  train- 
ing sessions  the  coniplete  nonlinear  [f(a,p)]  airframe  model  was  utilized  and 
departures  could  be  obtained.  For  the  final  data  runs  the  airframe  dynamics 
were  frozen  at  the  initial  a and  p values  to  provide  the  necessary  stationary 
control  process.  The  detailed  describing  function  and  performeuice  measures 
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are  contained  in  Ref.  31 • In  the  following  we  will  first  describe  results 
of  the  roll  tracking  task,  then  pitch,  and  finally  heading. 

A.  BOLL  PHOT  mASORBS  AID  MODEL 

The  open-loop  roll  attitude  dynamics  for  the  three  configurations  are 
shown  in  Fig.  46.  The  root  locus  and  Bode  for  the  uncoupled  Oo=  17.3°/Po  = 0° 
case  is  shown  at  the  top  of  Fig.  46  and  that  for  Configuration  1,  Oo=l8.8°/ 

Pq  = 3°,  is  shown  in  the  middle.  These  two  represent  the  nominal  A-7  aircraft 
dynamics  at  these  fli£^t  conditions.  The  bottom  Bode  and  root  locus  repre- 
sents the  Configuration  3,  ao=  l8.8°/Po  = 3°,  case  for  which  the  stability 
derivative  is  increased  from  -0.7  to  —1.6  and  Np  is  changed  fron  0 to  -K).3. 
The  latter  results  in  lateral  dynamic  pareuneters  which  closely  approximate 
those  of  the  uncoupled  ao  = 17.3°  case  but  with  the  lateral-longitudinal  coup- 
ling due  to  sideslip.  The  change  in  primarily  affects  the  longitudinal 
dynamics  as  will  be  shown  later. 

In  ■'•^e  following  we  shall  look  at  the  uncoupled  case  first  and  then 
exair*ne  the  effects  of  coupling. 

1.  ao  = 17.3°/3o=0® 

Figure  47  presents  the  open-loop  YpYg  sunplitude  (+)  and  phase  (□)  data 
points  calculated  from  the  describing  f mction  analyzer  (DFA)  readout  at  the 
five  forcing  function  frequencies.  The  plotted  points  are  the  average  of 
two  separate  runs  for  pilot  ML.  These  data  points  have  been  fit  with  a 
transfer  function  ^diich  provides  the  solid-line  amplitude  and  phase  shown. 

Also  shown  is  a reference  6 dB/octave  slope  line.  It  is  apparent  the  pilot 
has  adopted  equalization  so  that  YpYc  approximates  the  crossover  model  (Ref.  28). 
Crossover  is  obtained  at  I.5  rad/sec  and  35  deg  phase  margin. 

The  pilot  describing  function,  Eq.  6,  is  obtained  by  subtracting  the  known 
airframe  amplitude  and  phase  from  the  DFA  data  points  of  Fig.  47  and  curve 
fitting  the  resulting  points  (Fig.  48): 
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The  results  provide  an  excellent  curve  fit  and  also  validate  the  precision 
pilot  model  (Ref.  29).  However,  for  all  practical  purposes  the  low-freq\iency 
lag-lead  can  be  neglected  and  the  equalization  sinopllfled  to  a first-order 
lead  at  1.88  rad/sec. 

The  pilot  has  placed  the  lead  essentially  at  the  desired  crossover 
instead  of  at  the  roll  subsidence  lag  (O.55  rad/sec)  as  in  the  mojre  usual 
case.  In  this  instance  (Uq/cua  considerably  less  than  unity  and  provides 
lead  between  1 /Tr  and  the  desired  crossover.  This  reduces  the  need  for  lower 
frequency  pilot-generated  lead.  It  may  be  observed  from  Fig,  h6  that  the 
pilot  has  placed  the  series  lead  at  the  point  where  the  airframe  phase  lag  is 
approximately  155  <ieg  (i.e.,  at  the  "effective"  1/Tr). 

2.  Oo  = 10.0°/Po  = 3°#  Configuration  l 

The  measured  YpYc  data  points  and  the  subsequent  curve  fit  for  this  case 
is  shown  in  Fig.  49.  It  may  be  noted  by  CTmparing  Figs.  47  and  49  that 
the  YpYj.  for  this  case  is  essentieuLly  the  same  for  the  uncovq?led  case.  Sub- 
tracting the  controlled  element  amplitude  and  phase  fraa  the  DFA  measures  and 
curve  fitting  these  data  points  results  in  the  pilot  dynamics  shown  in  Fig.  50. 
The  transfer  function  for  this  data  fit  is  presented  as  Eq.  7; 


V = 5852(.66)(.76)  -.525s 

(-.05)[.5,  10](15) 

The  rig^t  half  plane  pole  at  “O.05  is  an  artifact  of  the  simplified  curve- 
fitting routine  for  the  increasing  amplitude  but  decreasing  phase  at  the 
lower  frequencies.  A more  accurate  phase  fit  would  be  obtained  at  low 
frequency  by  a model  containing  e''^/®j  however,  this  also  results  in  an 
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Figure  50. 


appreciable  cost  Increase  for  the  data  fit.  Beferrlng  hack  to  Fig.  k6,  it 
may  be  noted  that  the  airframe  has  considerably  less  low-frequency  gain  for 
this  configuratlcm  than  for  either  of  the  other  two  ccmflguratlons;  there- 
fore, the  pilot  adopts  a low-frequency  lag  to  Increase  the  YpYc  gain  at  the 
lover  frequencies.  A double  lead  is  then  required  in  the  vicinity  of  }/Tr 
to  obtain  adequate  phase  margin  (2^  deg)  at  the  crossover  of  1.6  rad/sec 
(Fig.  49). 

3 . a<,  = 18.0®/Po  = 3°»  Configuration  3 

The  describing  function  data  points  and  the  curve  fit  for  this  configuration 
are  shown  in  Fig.  ^1.  Again,  the  pilot  accomplishes  the  crossover  at  1.U  x ad/sec 
but  in  this  case  with  appraxlmately  UO  deg  phase  margin.  As  with  the  two  pre- 
vious cases,  the  pilot  adapts  equalization  to  achieve  a YpY(.  such  that  the 
crossover  model  applies.  Subtracting  the  known  Y^.  and  curve  fitting  the  result- 
ing amplitude  and  the  phase  data  points  provides  a pilot  model  given  by  Eq,  8 
and  shown  in  Fig.  52. 


In  this  Instance  zero  weight  was  assigned  to  the  lowest -frequency  phase  data 
point  (for  curve-fitting  routine).  Thus,  in  fitting  the  low-frequency  ampli- 
tude data  point  but  not  the  phase,  it  obtained  a left  half  plane  pole  at 
0.21  rad/sec.  As  might  be  expected  on  the  basis  of  similarity  of  airframe 
dynamics  between  this  case  and  the  uncoupled  case,  the  pilot  models  of 
Eqs.  6 and  8 are  essentially  the  same. 

1^.  Susnmry 

Figure  53  presents  a direct  comparison  of  the  YpYg  data  points  obtained 
for  the  three  configurations  and  Fig,  5U  presents  the  Yp  data  points.  It 
may  be  observed  that  the  data  exhibit  less  scatter  than  often  is  obtained 
for  multiple  runs  of  a single  configuration.  The  only  significant  difference 
between  the  configurations  is  in  the  low-frequency  region  where  it  appears 
the  coupled  lateral-longitudinal  airframe  dynamics  cause  the  pilot  to  use 
higher  gain  and  have  considerably  greater  phase  lag.  In  the  region  above 
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Figure  51*  Open  Loc^)  Ypq^c#  = l8.8°/Po  = Figure  52.  Yp_  Data  Fit,  Oo  = 18.8' 

Configuration  Pilot  ML  Configuration  5>  Pilot  ML 


0. ^  rad/sec  there  is  little  or  no  significant  difference  betveen  the  airfMme 
dynamics  en{d.pyed  and  little  reason  for  the  pilot  describing  function  measures 
to  differ  in  this  frequency  band.  The  conq^ling  had  little  Influence  on  the 
crossover  desired  or  achieved.  In  the  region  of  crossover  a single  pilot 
model  would  suffice  with  a first-order  series  lead  located  at  the  dutch  roll. 

The  major  effect  of  the  multiloop,  multi- input,  multi- output  task  (with  a roll 
command  forcing  function)  is  to  modify  the  pilot's  low-frequency  coupensation. 

B.  PIKE  FZXOI  IWOnfl  AID  ICOEL 

The  open-loop  pltdi  dynamics  for  the  three  controlled  element  configurations 
are  shown  in  Fig.  Again,  the  iqpeimost  Bode  and  root  locus  plots  represent 
the  uncoupled  case,  the  middle  plots  the  nominal  A-7  at  l8.8  deg  and  3 deg  Pq, 

and  the  lower  plots  the  Configuration  3 case  with  = 3 deg  and  = ->1 .6.  It 

may  be  noted  that  the  three  configurations  represent  three  different  values  of 
1/T0^  (from  0 to  -0.2).  There  is  also  considerable  difference  in  low-frequency 
amplitude  ratio  for  the  most  highly  coupled  configuration  (3)* 

In  the  following,  the  data  obtained  will  be  analyzed  on  two  bases.  First, 
we  will  neglect  the  contribution  of  the  roll  loop  closure  (see  Fig.  27)  and 
assume  the  controlled  element  is  described  by  the  single  open- loop  Yc  shown 
in  Fig.  33.  This  is  an  approximation  only  for  the  two  cases  having  lateral- 
longitudinal  coupling.  Second,  we  will  analyze  Configuration  1 assuming  the 
cooplete  multiloqp  structure  and  closing  the  inner  roll  loop  using  the  pilot 
model  obtained  in  Subsection  A- 2,  preceding. 

1 . Single-Loop  Modele 


a.  ao=17.?°/Pa.0° 


The  YpYg  data  were  fitted  for  one  set  of  measurements  for  Pilot  ML  and  are 
presented  in  Fig.  56.  These  are  identified  by  the  X and  □ in  Fig.  56.  These 
points  represent  an  average  of  two  runs  having  a maximum  of  ±2  dB  and  ±10  deg 
phase  from  the  data  point  plotted.  The  data  spread  for  two  runs  by  Pilot  JS  is 
indicated  by  the  ^ on  Fig.  56.  These  data  points  were  not  fitted j they  are 
presented  here  to  show  the  consistency  between  the  two  pilots.  Pilot  JS 


112 


ax 


Figure  56.  Open  Loop  Yp^Yc,  oq  = 17«5°/Po  = 0°  Figure  57-  Yp  Data  Point  Fit,  a©  = 17.3®/Po 


obviously  was  closing  the  loop  with  a lower  gain  and  therefore  had  a cross- 
over sometdiat  below  0,^  rad/ sec,  idiereas  Pilot  ML  achieved  a crossover  near 
2 rad/ sec.  A 6 dB  per  octave  reference  slope  line  Is  also  shown  on  Fig.  ^6. 

It  Is  apparent  that  both  pilots  closely  approximate  the  crossover  model. 

Pilot  JS  would  appear  to  be  placing  a first-order  lag  on  top  of  1 /T02  aod  a 
lead  near  a>op  as  might  be  expected  (see  Fig.  However,  the  transfer  func- 

tion fit  to  the  Pilot  ML  data  Indicate  he  Is  using  somewhat  more  complex  equali- 
zation. This  may  be  observed  by  subtracting  the  known  single-loop  (this  case 
Is  uncoupled)  Yq  anqplltude  and  j^ase  from  the  YpYc  data  points  of  Fig.  ^6  and 
fitting  them  as  shown  in  Fig.  57.  Uie  Yp  transfer  function  is  shown  in  Eq.  9: 


5ehi.%3)[.226.2.^]  -.258 

(-.08)(1.69)(10)[.5,  10) 


(9) 


Again,  this  conpares  well  with  the  precision  pilot  model  with  the  exception 
of  the  series  lead;  a low-frequency  lag-lead  is  enployed  to  Increase  the  low 
frequency  gain;  the  pilot  acts  as  a pure  gain  in  the  crossover  region;  and 
the  third-order  neuromuscular  system  lags  are  at  about  10  rad/sec.  The  main 
difference  is  the  more  complicated  series  equalization  which  consists  of  a 
second-order  lead  over  a first-order  lag.  From  the  amplitude  asymptotes  of 
Fig.  57  it  appears  that  this  complex  equalization  might  be  approximated  as  a 
.'dimple  first-order  lead;  however,  in  order  for  the  inflections  of  the  fitted 
curve  to  closely  match  the  actual  amplitude  ratio  and  phase  data  points  a 
lowly  damped  second-order  lead  is  required.  It  is  apparent  that  the  pilot 
is  precisely  inverting  the  short-period  characteristics.  The  current  preci- 
sion model  does  not  allow  for  such  a second-order  lead;  however,  the  precision 
model  was  also  obtained  with  Y^'s  having  much  more  ideal  dynamic  character- 
istics than  those  for  the  high  a case  employed  here.  That  is,  they  did  not 
have  second-order  dynamics  in  the  region  of  the  crossover  and  in  our  case 
the  pilot  is  achieving  crossover  right  at  the  short-period  mode.  Thus,  the 
series  equalization  might  be  expected  to  be  more  con5)lex. 

There  is  other  evidence  to  siq)port  the  existence  of  second-order  lead 
equalization.  A root  plot  of  the  early  Hall  data  (Ref.  10)  is  shown  in 
Fig.  58.  Two  data  points  at  roughly  1 .5  rad/sec  and  low  damping  ratio  are 
shown  to  lie  in  a region  in  which  double  lead  was  employed.  These  have 
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figure  58.  Con^rlson  of  Theoretical  Non- Equalized  Boundary  with  Hall  Measured  Cata 

(From  Ref.  10) 


•s- 


ratings  of  7 and  10.  A third  data  point  is  at  roughly  2.^  rad/sec  and  eLLso 
achieved  a rating  of  10.  It  may  be  noted  from  Fig.  55  that  the  short-period 
for  all  three  configurations  investigated  here  ■would  lie  in  aiqproximately 
the  center  of  the  double  lead  area  shown  in  Pig.  58.  The  Cooper -Harper 
ratings  of  8-^  obtained  for  the  three  configurations  in  this  simulation  are 
also  consistent  with  the  Hall  data. 

The  Shirley  data  (Ref.  11)  also  resul'ted  in  pilot  models  wi'th  second- 
order  lead.  The  example  transfer  functions  shown  in  Table  9 show  the  pilot 
is  capable  of  nearly  canceling  second-order  modes  of  the  controlled  element. 


TABLE  9 

E5CAMPLE  PILOT  TRANSFER  RJNCTIONS 
FR(»1  SHIRLEI  (Ref.  II) 


Yc 

Yp 

5 

1.16  e“'^®  (s  + .5)[s2  + 4s  + 5] 

[s^  + s + 5] 

(s  + .I3)(s  + 2) 

.0  = 

11  e~'5®  (s  + .4)rs2  + 3s  + 10] 

[s2  + 2s  + 10] 

(s  + .l)(s  + 2.32) 

ihe  precise  reason  for  the  Hall  euid  ^lirley  data  to  result  in  second-order 
leads  is  not  known.  For  the  task  employed  in  our  simulation  it  is  suspected 
that  the  pilot  learned  to  completely  suppress  the  short-period  dynamics  in 
order  to  pre'vent  departure.  That  is,  during  the  task  leeurning  period  with 
the  unfrozen  aerodynamics,  relatively  small  Increases  in  angle  of  attack  could 
result  in  depeu:ture.  Since  the  short-period  mode  is  predominantly  a and 
0 motion,  the  pilot  was  forced  to  control  a and  0 quite  tightly  to  suppress 
overshoots.  Such  complete  suppression  of  the  short-period  mode  would  then 
require  the  complex  equalization  shown  in  Fig.  57-  Howe'ver,  the  use  of  such 
equalization  would  not  be  suspected  from  the  open-loop  controlled  element 
dynamics  shown  in  Fig,  Ihere  are  also  no  known  pilot  rating  functionals 


for  secoM'Order  lead  equalization.  Urns,  It  would  be  difficult  Indeed  to 
predict  both  the  pilot  model  and  pilot  rating  for  pitch  attitude  control  at 
this  17.5  deg  a,  0 deg  p case. 

b . Ops  l8.8°/Po  = 3°i  Con^guration  1 

The  YpYc  amplitude  and  phase  data  points  for  this  case  are  shown  in  Fig.  39. 
Again,  the  X and  O represent  Pilot  ML  while  the  ^ represent  Pilot  JS.  Only 
the  data  points  for  Pilot  ML  have  been  fitted.  As  with  the  lateral  data,  the 
lowest-frequency  phase  i>oint  was  given  zero  weight  in  the  data  fit  for  economy 
of  conputation. 

For  this  case  the  crossover  model  seems  better  suited  to  the  lover  gain 
pilot  (JS).  Pilot  ML  adopted  a slope  somewhat  less  than  6 dB  per  octave. 

This  decreased  slope  has  been  obsei^d  previously  (Ref.  30 ) on  subcritical 
casks  where  the  controlled  element  has  a right  half  plane  pole.  Referring 
back  to  Fig.  55,  it  will  be  noted  that  this  configuration  does  have  a ri^t 
half  plane  zero  at  0.1  rad/sec.  Closing  the  0 &e  Jj^ives  the  open- 
loop  pole  from  the  origin  into  this  zero  and  therefore  results  in  a right 
half  plane  closed- loop  pole.  The  closed-loop  right  half  plane  pole  would 
normally  be  expected  to  result  in  a pitch  attitude  instability  that  would 
prevent  obtaining  describing  function  data.  However,  if  the  pitch  attitude 
loop  is  closed  tightly,  the  closed- loop  pole  is  essentially  driven  into  the 
open- loop  zero  and  there  is  very  little  pitch  modal  response.  Hie  small  pitch 
modal  response  together  with  the  low  frequency  provides  a very  long,  slow  pitch 
divergence  wh'  (ould  not  be  observable.  But,  the  instability  will  shov’  up 
in  another  mode.  x.ie  situation  is  analogous  to  tight  attitude  regulation  for 
precise  path  control  with  an  aircraft  on  the  backside  of  the  power  required 
curve.  The  i>ath  numerator  has  a right  half  plane  zero  which  gives  rise  to  an 
unstable  closed- loop  mode.  With  tight  attitude  regulation,  stable  path  control 
is  obtained  and  the  instability  shows  up  as  a speed  divergence.  For  our  high  a 
case,  the  right  half  plane  zero  arises  due  to  the  coupling  stability  derivative, 
N(J,.  From  the  previous  analysis  we  know  the  closed-loop  instabil5.ty  shows  up 
in  a lateral  mode  (nose  slice)  and  the  pilot  must  close  the  cp  6a  loop  to 
restabilize  it.  Therefore,  we  might  expect  the  slope  in  the  region  of  cross- 
over to  be  less  than  6 dB  per  octave.  The  end  result  is  that  the  pitch  attitude 
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control  task  is  slmlla~  to  the  uncotqpled  17.5°/0°  case  except  the  pilot  is 
r -w  required  to  close  a cp  5^  loop  to  decouple  the  latewil-longitudinal 
motion  eind  hence  is  working  harder. 

The  pilot  describing  function  is  again  obtained  by  subtracting  the  Yc 
from  the  YpY^  data  of  Fig.  59  and  curve  fitting  the  resulting  data  points. 

As  indicated  previously  we  will  first  assume  the  effective  Yq  is  adequately 
described  by  the  open  loop  0 5e  transfer  function  of  Fig.  55.  The  result- 
ing data  points  and  their  curve  fit  are  presented  in  Fig.  60.  The  pilot  model 
for  the  fit  is  given  in  Eq.  10: 


Y - 930.8(.2g)[.3 . 2.31  -.306s 

~ (.006)(5.i)(7.6)[.5,  10] 


(10) 


Again,  the  pilot  model  has  the  form  of  the  precision  model  with  a series 
equalization  consisting  of  a second-order  lead  over  a first-order  lag.  From 
Fig.  60  it  is  apparent  that  the  pilot  is  acting  as  a pure  gain  in  the  region 
just  below  crossover  and  is  applying  the  second-order  lead  to  cancel  the 
vehicle  short-period  mode.  In  this  case  the  effect  of  the  second-order  lead 
is  more  pronounced  because  of  the  greater  separation  between  it  and  the 
associated  lag  break  point. 

c.  g = l8.8°/p  = 3°t  Configiiration  3 

The  YpYc  data  points  euid  curve  fit  for  Pilot  ML  are  shown  in  Fig.  6l . 

Also  shown  are  the  raw  data  points  for  Pilot  JS.  Again,  the  slope  in  the 
region  of  the  crossover  is  less  than  that  which  would  be  obtained  with  the 
crossover  model  and  again  the  pitch  numerator  has  a right  half  plane  zero. 

In  this  case  it  is  at  0.2^i  rad/sec  and  provides  a significant  increase  in  the 
airframe  low-frequency  amplitude.  Also,  upon  closure  of  the  pitch  attitude 
loop,  the  unstable  root  is  driven  further  into  the  right  half  plane  and  thus 
the  mode  divergence  is  more  rapid.  The  unstable  mode  should  be  more  apparent 
to  the  pilots  and  hence  constrain  the  gain  and  equalization  necessary  to 
achieve  satisfactory  performance.  This  is  probably  the  reason  why  the  describ- 
ing function  data  points  for  both  pilots  are  in  closer  eigreement  fo:r  this  con- 
figuration than  for  the  previous  two. 


Figure  6l  • Open-Loop  Yp„Yc>  Qo  = l8.8°/Po  ~ Figure  62.  YpQ  Data  Fit,  Oo  = 18.8  /Po  ~ 

Configuration  5 Configuration  3>  Pilot  ML 


Again,  subtracting  the  c^n-loop  6 &e  airframe  anq^ilitude  and  phase 
from  the  YpY{.  data  points  of  Fig.  6l  gives  the  data  points  and  curve  dtown 
in  Fig.  62.  The  pilot  model  transfer  functicm  is  presented  in  Eq.  11: 


-.286s 

e 


(11) 


Here,  the  low-frequency  lag-lead  is  not  required  because  of  the  hi^er  gain 
of  the  Yq.  !Ihe  pilot  essentially  adopts  a pure  gain  out  to  the  short-period 
frequency  and  then  assumes  a second-order  lead  to  cancel  the  short-period 
resonance  x>eak. 

d.  Summary 

Figure  63  presents  a direct  comparison  of  the  YpY^  data  points  from  the 
three  configurations  for  Pilot  ML.  These  demonstrate  the  uniformity  the  pilot 
was  achieving  in  the  pitch  task  — coi;pled  versus  uncoupled. 

Figure  64  presents  a direct  comparison  of  the  Ypg  data  points  after  removing 
the  controlled  element  amplitude  and  phase.  As  with  the  roll  closure  these 
show  remarkable  similarity  in  the  region  of  crossover  and  above  and  could 
readily  be  fit  by  a single  model.  Again,  the  major  difference  is  in  the  low- 
frequency  band  with  the  greater  coupling  conditions  showing  the  lower  pilot 
anplitude  contribution. 

This  is  precisely  the  opposite  of  the  trend  shown  in  Fig.  54  where  the 
lateral  Yp^  for  the  coupled  case  showed  increased  low-frequency  gain.  It 
appears  that  the  pilot  is  achieving  a tradeoff  between  the  two  axes  — possibly 
to  effect  a decoipling  of  the  motions.  The  major  coupling  is  from  longitudinal 
to  lateral  via  Nq,.  therefore,  as  coupling  increases  (Nq,  increased),  the  pilot 
reduces  his  longitudinal  gain  emd  Increases  the  lateral  gain. 

In  all  cases,  whether  coupled  or  uncoupled,  the  pilot  is  closing  the  pitch 
loop  essentially  at  the  short  period  and  is  nearly  canceling  the  short-period 
peak  through  generation  of  a second-order  lead.  This  unusual  equalization 
has  apparently  been  adopted  throu^  the  necessity  to  avoid  excessive  angle  of 
attack  and  subsequent  nose  slice  divergence.  A clue  to  the  possible  tension 
of  the  pilot  in  pitch  control  is  the  ^aite  small  time  delay  (0.25  < t < 0.5  sec) 
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obtained  in  this  task  as  coo^iared  to  that  (0.5  < t < O.56  sec)  obtained  in 
roll.  This  nd^t  also  signal  that  the  pilot  closes  the  pitch  loc^  as  the 
inner  loop  and  roll  as  the  outer  locp. 

The  YpYc  anplitude  cmd  jAiase  did  not  fit  the  crossover  model  for  the 
two  coiq>led  cases.  Ihese  exhibited  considerably  flatter  amplitude  throu^^ 
the  crossover  region  and  are  consistent  with  previous  data  obtained  with  a 
controlled  elements  having  unstable  modes. 

2.  Multiloop  Model 

Because  the  coupling  from  longitudinal  motion  into  lateral  is  significant 
for  the  asynmetrlc  flight  cases.  Configuration  1 will  be  analyzed  by  multi- 
loop techniques.  The  extent  of  the  coupling  is  indicated  by  the  time  traces 
of  Fig.  65  which  are  taken  from  a pitch  command  data  run  with  Pilot  ML.  It 
will  be  noted  that  the  0 command  evokes  roughly  ±0.7  in.  of  longitudinal 
stick  response  and  a i>eak  iaa  of  about  ±2.9  deg.  This  then  couples  throu^ 
the  airframe  to  require  twice  as  much  lateral  stick  motion  and  a peak  beuik 
angle  of  over  ±10  deg.  On  the  other  hand,  the  coupling  from  lateral  motion 
into  longitudinal  motion  Is  relatively  insignificant.  Therefore,  we  will 
assume  the  cp  loop  is  the  inner  loop  and  determine  its  effect  on  the  outer,  9, 
loop.  The  inner  loop  is  assumed  to  be  closed  with  the  pilot  model  of  Eq.  7. 
The  effective  airframe  pitch  dynamics  are  then  obtained  as  shown  in  Eq.  12: 

® ^ A + M? 

Pqj  °a 

where 

= -2.91(-.i)(.21)(.595)[.279,  1.16] 

= 2.2[.25,  .57][.i8L,  1.92] 

Nbeba  = -1.29(.255)[.266,  .56] 

Y = ?852(.66)(.76)(-7.67)^ 

(-.03)[.5,  10](15)(7.67)2 

A = [.97J+,  .258][. 296,  1.05][.  185,  1.96] 

V.  - -2.9l(.2t*)(.60i*)(2.02)[.067  . .39][.12  . 1.87] 

~ (0)(  .U8)(1 .96) [.  18  , .U5][.  1 67  , 1 .7 J 1. 13^‘ , 2.07] 

12J+ 


The  Bode  asyaptote  for  the  resulting  effectlTe  pitch  airArame,  Yce,  is  shown 
in  Fig.  66.  Note  that  the  low-  to  mid-frequency  asymptote  closely  jreseables 
that  for  the  open-loop,  uncoupled,  17.3^0°  case  of  Fig.  33.  That  is,  the 
low-  and  mid-frequency  asyoqitote  of  Fig.  66  has  lowered  to  jresemble  that  of 
Fig.  33.  If  we  then  apply  the  pilot  model  derived  from  the  uncoupled  case 
(Eq.  9)  and  plot  Ype^ce  aoplltude  and  phase  at  the  forcing  function  frequen- 
cies, we  obtain  the  points  identified  by  X in  Fig.  67.  If  the  pilot  model 


Figure  66.  Bode  Asymptote  for  Pitch  Attitude  with 
Inner  Roll  Loop  Closed 


derived  from  the  coupled.  Configuration  1,  case  (Eq.  10)  is  applied,  we  obtain 
the  points  identified  by  O in  Fig.  67.  The  actual  measured  DFA  data  points 
for  the  multiloop  YpYc  are  shown  by  □ in  Fig.  67.  Note  that  while  all  points 
are  in  close  proximity  the  points  obtained  with  the  coupled  pilot  model  match 
the  measured  data  best  at  hl^  frequencies  while  the  points  obtained  with  the 
uncoupled  pilot  model  match  best  at  low  frequencies.  It  thus  appears  that 
pilot  closure  of  the  cp  -*■  6a  loop  effectively  decouples  the  lateral-longiutdinal 
dynamics  to  the  point  that  a pilot  model  derived  from  an  uncoupled  case  provides 
a somewhat  better  overall  fit  to  the  measured  multiloop  data.  However,  it  is 
also  apparent  from  Fig.  67  that  the  differences  are  minimal  and  therefore  the 
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Figure  67  • ^ce^p  Cotnparisons,  Pilot 
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single-loop  analysis  approach  of  the  preceding  Subsection  B-1  remains  valid. 
Further,  since  the  roll  closure  only  modifies  the  low-flrequency  pitch  dynamics 
it  may  again  be  concliided  that  the  closure  sequence  adopted  by  the  pilot  is 
pitch  as  inner  (ti^test)  loop  and  roll  as  outer  (lowest  gain)  loop. 

C.  HMDm  PILOT  NBflSllRES  AHD  MODEL 

Hie  heading  pilot  measures  were  obtained  by  inserting  the  forcing  function 
command  as  target  beuik  angle  but  with  the  target  fjree  to  displace  in  heading 
(i.e.,  = (g/UQs)9«p,  see  Fig.  27c).  Repeat  runs  were  made  recording  \^g  and 

9g  separately  to  facilitate  modeling  the  dynamics  of  each  loop.  The  heading 
and  roll  pilot  models  are  thus  subject  to  the  following  assumptions; 

1 . Hie  run- to-run  process  is  stationary. 

2.  The  roll  loop  is  closed  as  an  inner  loop  to  heading. 

3.  Coupling  from  lateral  to  longitudinal  motion  is  suffi- 
ciently slight  that  the  pitch  closure  does  not  alter 
the  lateral  dynamics  appreciably. 

The  data  points  obtained  from  the  T|rg  measures  for  five  runs  with 

the  Configuration  1 (l8.8°/3°)  airframe  dynamics  are  plotted  j.-:  Fig.  68.  It 
should  be  noted  from  the  legend  that  these  data  reflect  two  different  piloting 
techniques.  In  three  runs  the  pilots  were  instructed  to  use  only  ailerons  for 
lateral  control.  In  two  runs  they  were  allowed  to  use  rudder  if  they  desired. 
Rudder  was  not  considered  beneficial,  in  this  task  and  was  used  speuringly,  if  at 
all.  However,  this  appears  partially  responsible  for  the  considerably  greater 
scatter  than  obtained  in  the  previous  pitch  and  roll  measurements.  The  rather 
surprising  aspect  of  the  data  is  the  heading  crossover  above  1 rad/ sec.  The 
data  are  quite  consistent  in  this  factor.  Past  measures  of  heading  control 
in  landing  or  similar  tasks  have  generally  Indicated  crossovers  considerably 
below  1 rad/sec.  With  the  exception  of  three  amplitude  points  at  5 reid/sec 
the  data  is  in  good  agreement  with  the  crossover  model  (slope  6 dB/octave). 

The  effective  airframe  Yc^  of  Fig.  68  is  dependent  upon  the  inner  roll 
loop  closure.  As  previously  Indicated,  the  feedforw€u*d  cp^  path  of  Fig.  2T[c 
was  purposely  opened  to  determine  the  effect  of  this  path  on  overall  loop 
dynamics  as  measured  by  i^g.  The  data  points  for  Run  1102-18  of  Fig.  68 
reflect  the  absence  of  this  bank  angle  cue  €Uid  it  is  api>arent  that  the 
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removal  has  little  effect  on  the  heading  describing  function  data.  It  is 
therefore  concluded  that  the  pilot  pays  little  attention  to  target  bcuok  angle 
during  heading  tracking.  This  is  supported  further  by  cooparing  the 
data  points  obtained  from  the  q>g/q)  measures  for  the  three  loop  <p,  0) 
tracking  task  (Fig.  27c)  with  those  for  the  two  (cp,  9)  loop  task  (Fig.  27a). 
The  describing  function  data  for  these  two  cases,  two  runs  each,  are  plotted 
in  Fig.  69.  Ibe  three-loop  data  Indicate  considerably  less  an^litude  ratio 
at  the  lowest  frequency  and  much  less  phase  lag  at  the  two  lower  frequencies. 
This  is  consistent  with  the  pilot  not  tracking  in  roll  when  this  is  the  inner 
loop.  Therefore,  it  is  not  necesscury  to  adopt  the  low-frequency  lag-lead 
equalization  which  was  central  to  the  roll  tracking  task. 

Based  upon  the  comparison  of  data  in  Fig.  69,  a sinplified  inner-loop 
roll  pilot  model  was  selected  as  shown  in  Eq.  1?: 


Yp<p  = 


960( .6)  —.4s 

[.3,  10](15)  ® 


(13) 


This  reflects  the  first-order  lead  and  third-order  neuromuscular  lags  identi- 
fied in  Eq.  7 but  the  low-frequency  lag- lead  is  eliminated  and  the  gain  is 
reduced  by  roughly  a factor  of  four.  The  time  delay  has  also  been  reduced 
slightly  on  the  basis  that  it  would  be  necessary  for  the  pilot  to  reduce  his 
high-frequency  lags  if  the  outer  (\lr)  loop  is  to  be  stable  with  a closure  above 
1 rad/sec.  The  resulting  multiloop  structure  and  dynamics  are  shown  in  Fig.  70. 
The  open-loop,  amplitude  and  phase  plot  for  this  loop  structure  is  shown 

in  Fig.  71.  Also  shown  in  Fig.  71  are  the  describing  function  data  points  of 
Fig.  68  shifted  in  amplitude  to  correspond  to  an  outer  loop  Kp^  = O.625.  The 
derived  model  is  considered  to  provide  an  excellent  fit*  to  the  data  points 
measured  during  the  simulation. 

Again,  this  multiloop  model  reflects  a minimum  equalization,  minimum 
effort,  approach  on  the  part  of  the  pilot.  Apparently,  he  is  using  the  the 
inner  loop  primarily  to  damp  the  dutch  roll  nuisance  mode.  For  example,  the 
open-loop  is  O.296  whereas  the  closed  inner  loop  = O.58.  This  decreases 


*It  is  to  be  reemphasized  that  no  attempt  has  been  made  here  to  model  the 
low-frequency  phase  droop  due  to  the  "a  effect." 
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Figure  7O.  dynamic  and  Loop  Structure  for  Heading  Tracking  Task  Analysis 


the  amplitude  peak  at  dutch  zoll  sufficiently  that  the  outer  if  loop  can  be 
closed  tightly  (above  1 rad/ see)  without  the  necessity  for  additional  equali- 
zation. The  )|r  loop  with  its  hig^  dc  gain  then  provides  good  flight  path 
control . 

D.  SOMART  KSD  COKISJBIOKB 

Multiloop,  multi-input,  multi-output  pilot  meastirements  have  been  made 
for  an  air-to-air  tracking  task  with  aircraft  dynamics  representative  of 
in^iending  stall/departure.  These  measTores  have  allowed  identification  of 
the  effect  of  in^nding  stall/departure  on  pitch  control,  coupled  versus 
uncoupled  lateral-longitudinal  airframe  dynamics,  and  inner  versus  outer 
loop  applications  on  pilot  dynamic  parameters. 

The  crossover  model  was  found  to  apply  in  all  loop  closures,  including 
multiloop  situations,  except  when  the  controlled  element  exhibits  a right 
half  plane  pole.  In  the  latter  case  the  crossover  is  achieved  at  an  ampli- 
tude ratio  slope  somewhat  less  than  6 dB/octave. 

The  detailed  pilot  models  obtained  for  the  two-loop,  coupled  lateral- 
longitudinal  airframe  tracking  task  reflect  the  precision  pilot  model  form 
except  for  the  adoption  of  a second-order  series  lead  in  the  pitch  task. 

The  presence  (and  degree)  of  coupling  principally  influences  the  low-frequency 
lag-lead  equalization  adopted  by  the  pilot.  This  is  adjusted  within  and 
between  axes  to  effectively  decouple  the  airframe  motion  with  the  least  pilot 
effort  (equalization). 

Based  upon  the  open-loop  pitch  airframe  dynamics  used  here  it  would  not 
be  suspected  that  the  pilot  would  adopt  second-order  lead  equalization  for 
either  coupled  or  uncoupled  cases.  This  apparently  is  a result  of  the  require- 
ment to  prevent  pitch  (9  or  a)  overshoot  which  would  trigger  a departure.  Thus, 
the  pilot  adopts  equalization  based  on  crossover,  phase,  and  penalty  criteria. 
While  there  are  some  previous  pilot  parameter  measurement  data  which  also 
reflect  generation  of  second-order  lead,  these  are  insufficient  to  generate 
pilot  rating  functionals  or  to  even  predict  the  adoption  of  such  equalization. 

Considerably  more  data  must  be  obtained  before  a paper  pilot  type  model 
can  be  devised  for  the  pre-departure  flight  regime.  Most  past  effort  has  been 
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devoted  to  tasks  and  controlled  elements  where  the  airframe  has  relatively 
good  dynamics,  i.e.,  operation  well  within  the  nominal  flic^t  envelope, 
and  concern  is  related  to  identification  of  the  5-5  boundary.  This  high- 
angle-of-attack  task  has  involved  operation  at  (if  not  beyond)  the  i>er- 
missible  flight  envelope  boundary  with  attendant  airframe  dynamics  which 
have  been  shown  in  the  previous  sections  to  rate  in  the  6 to  9 region  at 
best. 

It  was  found  in  the  heading  tracking  task  that  the  pilot  was  not  making 
use  of  target  bank  angle  information  in  the  manner  usually  professed  by 
fighter  pilots,  i.e.,  match  target  bank.  One  possible  e3q>lanation  is  that 
this  technique  is  used  mainly  for  gross  maneuvers  and  that  once  in  a precise 
treu:king  position  the  pilot  switches  to  a "pointing”  mode  — which  this  simu- 
lation has  shown  can  accon^lish  a heading  crossover  of  greater  than  1 rad/sec 
even  with  relatively  poor  lateral  airframe  dynamics.  If  the  teurget  were  to 
make  a sudden  large  bank  angle  change  it  is  likely  the  pilot  would  immediately 
switch  to  a "match  bank  angle"  mode  to  minimize  either  error  or  reacquisition 
time. 
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SBCnOH  VI 


ccncLusicifs  AID  BECOMffinDAnars 


The  foregoing  analysis  and  simulation  has  traced  the  nose  slice  departure 
of  t e A-7  aircraft  throu^  a chain  of  events  Which  evolve  fron  unsymnetric 
fli^t  (directional  mistrim  or  miscoordination  in  maneuvers).  Biis  results 
in  sideslip  and  shed  vorticies  which  give  rise  to  aerodynamic  moments, 
and;fa>  ^ nonlinear  kinematic  coupling  also  occurs  between  sideslip  and 
yawing  or  rolling  of  the  aircraft  to  give  rise  to  pseudo-derivatives,  Zp 
and  Zj..  ihese  combined  coi^ling  effects  result  in  a right  half  plane  zero 
(1/T05  or  tue)  in  the  pi tch-attitvide- to- elevator  numerator.  Pitch  attitvide 
control  by  the  pilot  then  drives  the  aircraft  unstable.  The  instability 
starts  as  a first-order  divergence  in  coupled  lateral,  directional,  and 
longitudinal  motion.  However,  pilot  control  via  elevator  initially  restrains 
the  longitudinal  divergence  so  the  motion  appears  mainly  in  the  lateral- 
directional  degrees- of- freedom. 


This  somewhat  peculiar  control-divergence  phenomenon  is  analogous  to 
the  longitudinal  divergence  which  occurs  when  an  aircraft  is  on  the 
"backside”  of  the  thrust  required  curve.  In  the  latter  case  the  altitude 
numerator  has  a ri^t  half  plAne  (non-minimum  phase)  zero  and  control  of 
altitude  via  elevator  causes  a divergence.  However,  the  divergence  is  not 
reflected  in  altitude  (or  fli^t  path)  but  rather  in  speed. 


A quasi-linearized  five- degrees- of- freedom  model  containing  the  coupling 
terms  C Zp,  and  Zj.  and  with  only  six  lateral-directional  derivatives 
varied  as  a function  of  angle- of- attack  and  sideslip  provided  a good  simu- 
lation of  the  A-7  hi^  angle- of- attack  and  nose  slice  d^arture  characteristics. 


Based  on  the  nine  dynamic  configurations  examined,  it  is  concluded  that 
nose  slice  characteristics  and  pilox  comments  best  relate  to  the  longitu- 
dinal parameters  l/^e^  “0*  These  are  strongly  influenced  by  the 

derivative  N^.  The  conventional  lateral  parameters,  and  certainly 

contribute  to,  but  are  not  the  key  factors  in,  the  A-7  nose  slice.  However, 
it  appears  that  when  these  lateral  parameters  become  negative  at  the  same 
time  the  longitudinal  parameters  are  in  the  ri^t  half  plane  the  rapidity 
and/or  severity  of  nose  slice  d^arture  is  increased. 
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In  the  fixed  base  simulation  the  nose  slice  onset  could  only  be  detected 
vhen  the  target  aircz«ft  was  in  view.  Pilot  attempts  to  oppose  nose  slice 
lead  to  saturation  of  lateral>directional  control.  If  system  dynamics  or 
control  chau:acteristics  do  not  change  appreciably  as  departure  is  approached, 
the  lack  of  cues  can  Induce  the  pilot  to  continue  control  activity  until 
the  departure  is  well  developed  and  thus  lead  to  a more  severe  dqE>arture. 

The  configuration  having  the  best  lateral-directional  characteristics  but 
also  the  least  changes  in  c^namics  with  a was  downgraded  severely  because 
of  inadequate  warning.  In  this  respect  a steady  (graceful)  degradation  in 
handling  and  a slow  but  firm  nose  slice  provided  the  best  warning  and  time 
to  react. 

The  principal  means  of  recovery  is  to  open  the  pitch  attitude  loop,  i.e., 
let  go  of  the  stick,  at  the  first  indication  of  nose  slice.  The  use  of 
lateral-directional  control  generally  aggravates,  rather  than  assists, 
recovery.  One  key  to  successful  recovery  with  minimal  pilot  workload  or 
conpensatiai  is  that  the  aircraft  be  trimmed,  stick  force  zero,  at  an  an^e- 
of-attack  for  which  the  lateral-phugoid  and  dutch  roll  are  stable.  Therefore, 
the  open-loop  parameter  is  a key  factor  in  recovery. 

Data  and  pilot  ratings  obtained  for  the  stick  release  only  recoveries 
probably  relate  mostly  to  1 g,  strai^t-ahead  stall  maneuvers.  With  few 
exceptions  these  involved  minimum  penetration  into  stall/departure,  were 
readily  recoverable  control  free,  and  received  Cooper-Harper  ratings  of 
5 to  6 from  one  pilot  and  5 to  7 for  the  other  (using  stick  push).  Data 
and  ratings  obtained  for  all  other  types  of  recoveries  (e.g.,  B,  C,  D) 
probably  relate  better  to  actual  tracking  and  maneuvering  situations  in 
which  inadvertent  departure  — or  delay  of  departure  — is  encountered. 

These  cases  consistently  resulted  in  Cooper-Harper  ratings  of  6 to  10. 

The  very  significant  difference  between  these  two  sets  of  ratings  may 
reflect  actual  fli^t  situations  in  which  stall  and/or  spin  characteristics 
obtained  from  intentional  entries  (including  application  of  pro-spin 
controls)  tend  to  be  considerably  milder  than  those  obtained  from  inadvertent 
entries. 

The  resulting  pilot  ratings  could  not  be  correlated  with  any  single 
flying  quality  parameter.  The  departure  characteristics  were  most  dependent 
upon  closed  loop  (i.e.,  pitch  numerator)  parameters  with  some  influence 
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frcm  open  loop  paraaeters.  Where  recovery  vas  possible  it  vas  aost  dependent 
vpon  alrfrcuse  open  loop  parameters  (l>e.,  denominator  and  stick  trim)  but 
the  possibility  of  successful  recovery  vas  also  hl^ily  dependent  upon  pilot 
tecdmique,  agressiveness^  sensitivity  to  cues,  etc. 

All  configurations  investigated  here  exhibited  good  pitch  control  and 
response  characteristics.  If  a configuration  exhibiting  plt<di-up,  loss 
of  longitudinal  control  effectiveness,  or  ai^  other  impairment  of  longl> 
tudinal  control,  were  to  be  enployed  the  results  mi^t  be  considerably 
different  from  those  obtained  here. 

The  resxilts  of  this  study  indicated  the  conventional  closed-loop  pilot/ 
vehicle  analysis  approach,  while  appropriate  for  pre-departure,  is  not 
appropriate  for  acttial  depsirture  and  recovery.  Immediately  prior  to  depar- 
ture the  lateral-directional  control  (if  employed  at  all)  tended  toward 
steady  or  trim  type  application  against  the  motion  of  the  nose  slice. 
Following  departure  the  preferred  method  of  recovery  was  to  release  all 
cmtrols.  ihus,  the  pilot  provided  no  inputs  and  the  recovery  was  achieved 
through  open-loop,  airframe- alone  stability  and  daaping.  Again,  if  control 
was  applied  during  the  recovery,  the  preferred  technique  was  to  hold 
constant  forward  stick  to  assure  decrease  in  a and  to  prevent  pitch  up 
due  to  inertia  coupling.  As  the  aircraft  motion  subsided  the  stick  was 
returned  to  neutral.  Thus  a closed  loop  pilot  model,  per  se,  does  not 
appear  appropriate  for  departure/recovery  application. 

In  the  pre-departure  phase  where  closed-loop  considerations  are  valid, 
it  was  found  that  the  pilot  adopts  a second  order  lead  in  pitch  attitude 
control  — apparently  in  an  attenpt  to  cancel  short  period  pitch  attitude 
excursions  which  mi^t  trigger  departure.  This  second- order  complex  leeid 
generation  has  been  observed  before;  however,  there  are  insufficient  pilot 
rating  data  to  establish  a functional.  The  vast  amount  of  data  obtained 
for  first  order  lead  generation  is  not  applicable.  Much  more  data  is 
required  from  departure  type  conditions  before  trends  can  be  established 
between  pilot  parameters,  pilot  ratings,  and  performance  measures. 

A key  proVxem  with  ratings  exists  because  loss- of- control  is  involved. 
That  is,  acconplishment  of  the  primary  task,  \diatever  it  may  be,  has 
resulted  in  the  aii  craft  approaching  stall/ departure . If  the  task  is 
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pxirsued,  it  presumably  will  result  in  dqtarture  — loss  of  ccsitraiL. 
Departure  prevention  then  requires  abandoment  of  the  primary  task  and 
shifting  full  priority  to  safety  of  flight.  Thus,  Level  3 flying  quali- 
ties are  involved  — and,  in  fact,  the  question  is  establishment  of  the 
Level  3 boundary,  i.e.,  loss  of  control.  On  the  Cooper-Harper  scale  this 
is  the  boundary  between  9 10  (there  is  no  1^  aJi  airplane  will 

depart,  then  there  is  no  question  regarding  the  existence  of  the  Cooper- 
Harper  boundary,  but  merely  a questi<xi  of  how  much  of  the  potential  fli^t 
envelope  must  be  traded  for  fli^^t  safety.  !Dms  a mission  performance, 
usable  load  factor,  flig^it  safety,  or  some  related  index  becomes  paramount. 

In  a limited  investigation  of  ving-rock,  an  attempt  to  change  the 
d^arttire  characteristic  of  the  A-7  aircraft  from  nose  slice  to  wing  rock 
proved  partially  sucessful.  Scaling  of  five  key  stability  derivatives 
(jf*,  .fp,  N^,  Np,  and  to  resemble  those  of  the  F-4  did  result  in  an 
oscillatory  characteristic;  however,  it  did  not  have  the  predominate  roll 
characteristic  of  wing  rock. 

This  investigation  has  only  been  a first  cut  at  a most  ccnplex  problem 
area.  However,  it  has  yielded  considerable  insist  to  the  nose  slice 
phenomenon.  It  is  recommended  that  the  approach  be  continued  and  expanded 
in  a number  of  ways: 

1.  Investigate  a larger  population  of  aircraft  to 
validate  and/or  refine  the  key  parameters  identi- 
fied  here. 

2.  Develop  literal  expression  for  the  coupled  lateral- 
longj.tudinal  transfer  functioi  parameters  which 
result  from  unsymmetric  fli^t. 

3.  Investigate  possible  limits  of  (or  for)  rudder 
maneuvering  at  hi^  angle-of-attack. 

4.  Investigate  the  influence  of  single-  and  milti- 
axis,  hi^-gain,  hi^-authority  augmentation 
systems. 

5.  Expand  the  pilot  describing  function  data  base 
for  these  coupled,  limit- of- ccaitrol  situaticais. 

6.  Expand  the  simulation  to  six  degrees- of- freedom 
and  investigate  the  influence  of  additional 
aerodynamic  nonlinearities  (e.g.,  0[^,  Cjn^,  etc.). 
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It  would  also  be  hl^ly  desirable  to  investigate  the  Influence  of  motion 
cues  on  pilot  detection  of  departure  onset  and  assessment  of  departure- 
recovery  dynamic  characteristics.  But  unfortunately  the  raqpld  large  yaw 
and  lateral  acceleration  of  nose  slice  may  not  be  compatible  with  practical 
displacement  limits  of  most  moving  base  simulators.  The  motion  wai^out 
necessary  to  prevent  hitting  displacement  limits  mi^t  negate  completely 
the  desired  acceleration  cue. 
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